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The  purpose  of  this  study  was  to  investigate  the  interaction  of 
plane  elastic  waves  with  a  thin^  hollow,  cylindrical  shell  embedded  in  an 
elastic  medium. 

The  cylindrical  shell  is  considered  to  be  elastic,  isotropic, 
homogeneous,  and  of  infinite  length.   It  is  surrounded  by  an  elastic,  iso- 
tropic, and  homogeneous  medium  whose  motions  conform  to  the  ordinary  theory 
of  elasticity.  A  plane  stress  wave,  either  dilatational  or  shear  with  a 
step  variation  in  time,  whose  wave  front  travels  in  a  direction  perpendic- 
ular to  the  cylinder  axis,  envelops  the  shell.   Later,  a  Duhamel  integral 
is  used  to  study  other  wave  shapes  for  the  incident  stress- 

The  response  of  the  shell  is  studied  by  expressing  the  two  com- 
ponents of  displacement,  radial  and  tangential,  in  terms  of  Fourier  series, 
each  term  of  which  is  called  a  mode.   The  equations  of  motion  of  the  shell 
in  vacuo  are  derived  from  expressions  giving  the  strain  and  kinetic  ener- 
gies due  to  generalized  external  forces .   Forces  on  the  shell  result  from 
the  stresses  in  the  medium  at  the  boundary „   Stresses  in  the  medium  are 
taken  to  be  the  sum  of  the  stresses  due  to  the  incoming  stress  wave  expressed 
in  terms  of  Fourier  series  whose  coefficients  are  known,  and  those  due  to 
the  reflected  and  diffracted  effects  expressed  In  terms  of  a  pair  of  dis- 
placement potentials  representing  waves  diverging  from  the  axis  of  the 
shell. 

The  equations  to  be  solved  consist  of  two  pairs  of  coupled  integro- 
differential  equations  in  the  generalized  coordinates  of  the  shell  and  the 


displacement  j:.tentials.   By  use  of  a  digital  computer  they  are  solved  mode- 
wise  by  a  step-by-step  iterative  integration  technique  known  as  the  Newmark 
Beta  Method,  with  which  values  for  the  potential  functions^  and  the  accelera- 
tions, velocities,  and  displacements  of  the  shell  are  determined.   Stresses 
in  the  shell  are  found  from  the  displacements,  and  the  values  of  the  potential 
functions  permit  determination  of  stresses  for  any  point  in  the  medium. 

Although  the  equations  are  written  to  include  an  infinite  number 
of  modes,  only  the  first  three  modes  are  considered  in  detail.   Tne  computed 
solution  is  compared  with  values  obtained  from  a  series  expansion  of  the 
equations,  which  is  valid  for  short  times,  and  with  the  static  solution  based 
on  the  theory  of  elasticity  to  which  the  general  solution  should  approach 
asymptotically.   In  addition,  the  results  of  two  particular  problems  are 
compared  with  results  given  in  another  study. 

Numerical  solutions  are  obtained  to  determine  the  effect  of  the 
several  parameters  which  describe  the  relative  physical  properties  of  the 
shell  and  medium.   Results  are  presented  in  tabular  and  graphical  form. 
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CHAPTER  I 
INTRODUCTION 

1.1  General  Remarks 

The  problem  of  designing  underground  protective  structures  to 
resist  the  effects  of  nuclear  weapons  has  "become  increasingly  important  in 
recent  years  with  the  development  of  modern  weapons  whose  destructive  capacity 
is  overwhelming.   Engineers  in  this  field  are  hampered  to  a  great  extent  by  a 
lack  of  theoretical  information  on  how  structures  in  media  such  as  soil  or 
rock  behave  when  subjected  to  dynamic  loads.   Even  for  static  loads  alone, 
much  of  the  design  practice  today  is  of  a  semi-empirical  nature. 

When  a  nuclear  explosion  occurs,  stress  waves  are  transmitted 
through  the  air  and  ground.   How  are  they  transmitted  and  how  are  they  modified 
by  the  presence  of  a  structure  embedded  in  the  medium?   How  does  the  structure 
respond? 

The  purpose  of  this  report  is  to  study  one  aspect  of  the  problem, 
the  interaction  between  the  medium  and  stmcture. 

1.2  Statement  of  Problem 

The  problem  considered  here  consists  of  analyzing  the  elastic 
response  of  a  hollow  cylindrical  shell  (tunnel  lining)  embedded  in  an  elastic 
medium  when  subjected  to  an  incident  plane  stress  wave  traveling  in  a  direction 
perpendicular  to  the  axis  of  the  shell.   Some  questions  with  which  this  problem 
may  be  associated  are:   Do  tunnel  linings  in  contact  with  rock  afford  a  measure 
of  protection  significantly  higher  than  an  unlined  opening?   What  magnitude 
and  time  variation  of  displacement,  velocity,  and  acceleration  would  equipment 
mounted  within  such  a  structure  be  subjected  to?   How  are  stress  waves  within 


the  medium  modified  in  the  vicinity  of  the  shell?   This  study  was  conducted 
in  an  attempt  to  find  some  qualitative  and  quantitative  answers  to  these 
questions  within  the  limitations  imposed  by  the  assumptions  noted  below. 

1.3  Basic  Assumptions 

The  cylindrical  shell  is  considered  to  be  of  infinite  length,  and 
is  embedded  in  an  elastic  medium  of  infinite  extent  in  all  directions.   A 
plane  stress  wave  whose  front  travels  in  a  direction  perpendicular  to  the 
cylinder  axis  envelops  the  shell.   Strains  parallel  to  the  axis  in  both  the 
mediiim  and  shell  are  assumed  to  vanish;  thus,  since  each  cross  section  of  the 
shell  is  exactly  similar  to  every  other,  the  problem  is  reduced  to  one  of 
plane  strain. 

In  the  mathematical  development  of  the  problem  certain  basic 
assumptions  were  made,  the  most  important  of  which  are  given  here  with  a  few 
explanatory  remarks : 

(1)  The  medium  is  considered  to  be  homogeneous,  isotropic,  and 
linearly  elastic.   This  implies  that  the  ordinary  theory  of  stress  wave 
propagation  applies.   In  view  of  the  non -homogeneous,  non-isotropic,  and 
non-elastic  characteristics  of  most  materials  encountered  in  nature,  this 
is  a  severe  limitation;  however,  current  theories  of  stress  propagation 
through  such  media  have  not  advanced  to  the  stage  where  this  limitation  can 
be  readily  overcome.   In  the  case  of  some  rocks,  though,  this  assumption 
may  be  justified. 

(2)  The  material  in  the  shell  is  also  considered  to  be  homogeneous, 
isotropic,  and  linearly  elastic.   Generally  speaking,  this  assumption  is  valid 
for  values  of  stress  below  the  so-called  proportional  limit  of  materials 
commonly  used.   In  addition,  the  thickness  of  the  shell  relative  to  its  radius 


Is  assumed  small;  this  permits  expression  of  all  stress  components  of  the 
shell  in  terras  of  a  function  which  describes  the  deflection  of  its  middle 
surface.   This  deflection  must  satisfy  a  linear  partial  differential  equation 
with  the  appropriate  boundary  conditions. 

(3)   The  incident  stress  wave  considered  is  either  a  plane  dilata- 
tional  or  a  plane  distortlonal  (shear)  wave.   Under  actual  conditions,  both 
waves  are  propagated  with  the  shear  lagging  the  dilatational  wave.   The 
combined  effect  for  elastic  conditions  may  be  determined  through  the  principle 
of  superposition. 

(k)     The  radial  and  tangential  particle  velocities  of  the  medium 
at  the  boundary  are  equal  to  that  of  the  shell.   This  is  the  continuity 
relation  insuring  that  the  shell  and  the  medium  are  in  contact  with  no 
relative  slip  occurring  at  the  boundary. 

(5)  Any  additional  mass  within  the  shell  is  assumed  to  be  dis- 
tributed symmetrically  about  the  axis.   The  significance  of  this  assumption 
is  found  in  the  development  of  the  equations  of  motion  to  account  for  any 
additional  mass  located  within  the  shell. 

Other  assumptions  a  re  presented  in  the  formal  development  of  the 
mathematical  expressions  used  to  describe  the  behavior  of  shell  and  medium. 

1.4  Method  of  Approach 

The  two  components  of  shell  displacement,  radial  and  tangential, 
are  written  in  terms  of  Fourier  series  from  which  expressions  giving  the 
strain  energy  and  kinetic  energy  of  the  shell  in  vacuo  are  derived.   The 
equations  of  motion  are  derived  from  Lagrange's  equations  in  terms  of  the 
displacQTient  functions  and  forces  acting  on  the  shell.   Forces  on  the  shell 
result  from  the  stresses  in  the  medium  at  the  boundary.   Stresses  in  the 


mediiom  are  taken  to  be  the  sum  of  the  stresses  due  to  the  incoming  stress  wave 
expressed  in  terms  of  Fourier  series  whose  coefficients  are  known,  and  those 
due  to  the  reflected  and  diffracted  effects  expressed  in  terms  of  a  pair  of 
displacement  potentials  representing  waves  diverging  from  the  axis  of  the 
shell.   The  form  of  these  potentials  as  derived  by  Lamb  [k)    ,    is  in  terms  of 
sine  and  cosine  series^  the  nature  and  treatment  of  which  has  been  studied 
by  Paul  (6)0 

The  equations  of  motion,  described  as  a  pair  of  coupled  integro- 
differential  equations,  are  solved  modewise  using  a  numerical  technique 
known  as  the  Newmark  Beta  Method  (5)  which  permits  determination  of  the 
coefficients  of  the  potentials  and  values  of  acceleration,  velocity,  and 
displacement  of  the  shell. 

The  solution  obtained  is  compared  with  A^alues  obtained  from  a  series 
expansion  of  the  equations,  which  is  valid  for  short  times,  and  with  the  static 
solution  based  on  the  theory  of  elasticity  which  the  machine  solution  should 
approach  asymptotically. 

1.5  Previous  Work 

The  problem  stated  above  has  been  the  subject  of  a  recent  report 
by  Baron  (l),  whose  analysis  consists  of  first  solving  for  the  displacements 
caused  at  the  boundary  of  an  unlined  cylindrical  cavity  subjected  to  a  plane 
stress  wave.   This  is  done  through  an  integral  transform  approach,  the  solu- 
tion of  the  transformed  equations  being  expressed  using  Hankel  functions. 
The  evaluation  of  the  inverse  transform  is  accomplished  only  with  great 
computational  effort.   Values  of  displacements  obtained  are  then  used  as 
influence  coefficients  in  determining  the  displacements  of  the  shell. 


*  Numbers  in  parenthesis  refer  to  the  corresponding  entry  in  the  Bibliography, 


Solutions  olDtained  for  two  different  sets  of  parameters  by  Baron  are  compared 
in  Section  ^4-0  6  with  results  obtained  by  the  method  of  solution  outlined  in 
this  report o 

The  study  by  Paul  (6)  consisted  of  analyzing  the  effect  of  a  plane 
stress  wave  incident  on  an  unlined  cylindrical  cavity  In  an  elastic  medium. 
The  reflected  and  diffracted  waves  are  described  in  terms  of  displacement 
potentials  which  represent  outgoing  shear  and  dilatational  waves o   A  method 
was  developed  for  determining  values  of  these  potentials^  and  a  similar 
method  is  used  in  this  report. 

I06  Notation 

Notation  is  defined  throughout  the  text  where  it  first  appears; 
however,  the  following  list  sumjnarizes  the  main  uses  of  certain  symbols.   In 
discussions  of  special  topics  other  meanings  may  be  ascribed  to  the  symbols, 
at  which  time  they  will  be  redefined. 

A  =  Area  of  cross  section  of  the  shell,  per  unit  length. 

A  ,  B  ,  C   -     Coefficients  of  Fourier  series  for  stresses  in  the 
n'   n^   n 

medium. 

A    ,    B   ,    C        =     Coefficients   of  Foijrier  series   for  total  stresses 
n       n       n 

in  the  medium . 

AM  ,  BM     =  Weighting  factors, 
m^   m         too 

a  5  b       =  Generalized  coordinates  for  the  displacements  of 
n   n 

the  medium  at  the  boundary.. 

a   ,  b      ■■=  Generalized  coordinates  for  displacements  of  the 
ns   ns        ,  -, , 
shell.. 

c, ,  c       =  Velocities  of  wave  propagation,  dilatational  and 
shear,  respectively. 

d  =  Distance  from  neutral  axis  of  shell  to  its  outermost 

fiber. 


E,  E        =  Moduli  of  elasticity  for  the  medium  and  shell, 

respectively o   Bar  over  the  symbol  refers  to  the 
plane  strain  modulus, 

e  =  Volumetric  strain. 

F,  G        =  Generalized  coordinates  of  the  displacement 

potentials. 

1  =  Moment  of  inertia  of  the  shell,  per  unit  length. 

k  =  Parameter  which  relates  to  the  shape  of  a  time 

dependent  stress  wave. 

k  =  R/r,  ratio  of  radii. 

r  '  ' 

k  =  q.^/q.,,    ratio  of  velocities. 

c  2'  1 

m  =  Mass  of  the  shell  per  unit  surface  area, 

m'  =  Additional  mass  within  the  shell, 

n  =  Mode  number » 

=  Generalized  forces. 

=  Radius  of  the  shell.. 

=  Multip.lying  factors » 

r,  9        =  Polar  coordinates „ 

,T  =  .Kinetic  energy  of  the  shell. 

T'  =  Kinetic  energy  of  additional  mass  within  the  shell. 

t  =  Thickness  of  the  shell;  also,  time. 

U  =  Strain  energy  of  the  shell. 

u  =  Displacement  vector. 

u,  V        =  Components  of  the  displacement  vector  in  the  radial 
and  tangential  directions,  respectively. 

u  ,  V       =  Radial  and  tangential  displacement  components  of 
the  shell. 

u  ,  u       =  Components  of  the  displacement  vector  in  the  x  and  y 
directions,  respectively. 

X,  y,  z     -     Rectangular  coordinates. 


^n^ 

\ 

R 

R    , 

m' 

% 

a   E 
sn  s 
Q!  =  — =r J  non-dimensionalized  radial  component  of 

p    displacement. 

Td   E 
sn   s 
P  =  — ,    non-dimensionalized  tangential  component 

p    of  displacement. 

e  ,  e       =  Strains  along  the  x  and  y  axes,  respectively. 
X   y 

e„         =  Circumferential  strain, 
o 

t,  =  A  variable  of  integration. 

E 
T]^         =  g^  ,  a  parameter. 

T].  =  —  ,  a  parameter. 

T)  =  ~  ,    a  parameter. 

■n  =  —  _,  a  parameter. 

t  K 

Ps 

Ti  =  —  ,  a  parameter. 

P  P 

•n  =  -^-^ — /■:;  ^    \     ■       ,   a  parameter. 

Q,   r  =  Polar  coordinates. 

Q^  =  Position  angle  of  the  incident  wave. 

k  =  Curvature  of  the  shell. 

X,   [i  =  Lame'  constants. 

V  =  Poisson's  ratio. 


V 


V 

1-v 


I  ;  ^       =  Arg\;uiients  of  the  functions  F  and  G^  respectively. 

p  =  Mass  density  of  the  medium. 

p  =  Mass  density  of  the  shell. 

0,0  =     Normal  components  of  stress  parallel  to  the 

X  and  y  axis^  respectively. 

a  =     Shearing  stress  component. 


a   ,  0^^    =  Radial  and  tangential  normal  stresses  in  polar 
rr   yy 

coordinates.  Additional  subscripts  s  and  m, 

when  used,  refer  to  the  shell  and  medium. 

a  ^         =  Shearing  stress  in  polar  coordinates . 

a  ,         =  Bending  stress  in  the  shell, 
sb 

a  =  Maximum  stress  in  the  shell, 

smax 

a  =  Amplitude  of  incident  dilatational  wave. 

a  =  Amplitude  of  incident  shear  wave. 


=  Non-dimensionalized  time  in  the  case  of  the 
incident  dilatational  wave. 


"^^        -     Non-dimensionalized  time  in  the  case  of  the 
incident  shear  wave. 

Cp,  \1;        =  Displacement  potentials,  dilatational  and  shear, 
respectively . 


CHA.PTER  IT. 
MSIC  EQUATIONS 

2.1  Equations  for  the  Medium 

The  differential  equation  of  motion  of  a  particle  in  a  homogeneous^ 
isotropic,  and  linearly-elastic  medium  in  terms  of  its  displacement  vector  u 
is  given  by  Kolsky  (3)  in  the  form 

(\  +  2ii)  V  V  •  u  -  ^lV  X  (V  X  u)  =  p  ^  (2-1) 

^t^ 

where  X   and  \i   are  the  Lame  constants,  and  p  is  the  density  of  the  medium.   The 

vector  u  may  he  expressed  as  the  sum  of  two  displacements,  the  fi\v-OLci«e,\At   of 

a  scalar  potential  and  the  curl  of  a  vector  potential 

u  =  Vcp+Vxf  (2-2) 

Here  cp  is  a  potential  giving  rise  to  an  irrotational  displacement  and  Y  a 
potential  leading  to  an  equivoluminal  displacement. 

If  the  wave  equations  (2-3)  are  satisfied  and  u  expressed  as  in 
(2-2),  the  equations  of  motion  (2-l)  are  automatically  satisfied. 

2     1  ^ 

c   dt 


(2-3) 


2      1  S^ 

2  ^,2 
c  dt 


where  the  velocities  of  wave  propagation  are 


\+2ll 

^1  ^  y  p 

(2-4) 


2    Vp 

for  the  dilatational  and  shear  waves,  respectively. 


10 


Since  the  problem  is  essentially  a  two-dimensional  one  wherein  only  the 
rotation  about  the  cylinder  axis  is  considered,  the  expression  for  the 
propagation  of  shear  waves  can  he  written  in  terms  of  a  scalar  potential 
function 

V^  4r  =  A  ^   ;    with  ^  =  Y^  (2-3a) 

cp  and  ^   are  functions  of  x,  y,  and  to   The  components  of  the  displacement 
vector  u  can  be  expressed  as 

u  3.4i  (2-5) 


X   ox 


"¥ 


y  ~  ^   ^x 


Stress  components  are 

a       =  Xe  +   2u,e 

XX  x 

a   =  |J,Y  (2-6) 

xy   "^  xy  ^   ' 

a   =  Xe  +   2Lie 

yy       y 

where  e  represents  volumetric  strain,  e  and  e  strain  along  the  x  and  y  axis 

respectively,  and  y   the  shearing  strain.   The  strain  components  for  small 

xy 

displacement  are 


e 

X 

au 

X 

^x 

e 

y 

au 

7 

^xy  = 

du          ^u 
qy          ox 

(2-6a) 


In  terms  of  the  potential  functions,  the  equations  for  the  stress  components 
become 
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XX 


xy 


yy 


2  2  2  2 

^     /^cp        ^CDs         ^      /d  cp  ^_^_J[n 


bx         by 


\  2     ax^' 
ox 


=  1^ 


:^2  :s2 ,         ^2 , 


:.2  :,2  .2  ^2 , 

M^  +  ^)  .  2,(^  -  IJ-) 


(2-7) 


In  polar  coordinates 
.2 


2  ^  £_   1  ^   J:.  £_ 

'    .  2  ^  r  ar  "^  2  .„2 
or  r  oc^ 


and  the  wave  equations  expressed  in  polar  coordinates  are 


.  2   r  ^r    2  .„2 
or  r  dd 


2  ;^^2 
C-,  ot 


.  2  "*"  r  "SF  "^  2  -,.2  ~  2  -,,2 
or  r  o0    c   ot 


(2-8) 


Displacements  in  polar  coordinates  become 


u  - 


)r   r  ^0 


V  = 


1  ^   b^ 
r  ae  "  ^r 


(2-9) 


where  u  and  v  represent  radial  and  tangential  components  of  displacement, 
respectively- 


Stress  components  may  be  written  as 

'^   1  b^^       1_  bi 

^  2  ^  r  5r5e  ~     2  ^ 
-or  r    -• 


^^^  =  (^  +  2^) 


\ 


2,  -! 


1  a  cp    1  aijf    ^  jf 

2+  r  '§9  ■  ar^_ 


be 


^e  =  t^ 


2  ^%  2_  ^   ^   _^  ^   1  ^ 

r  arde  ~2b9'^2^'^        o  +  ^:^„ 
r       or 
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/  X+2u,  >. 


^   1  a  cp   1  b^        b  ^ 
br   ^  r\Q2  "^  r  '5^  "  ^r^Q 


r^  ^9^   r  ^rj 
+  \ 


b%       1  a^^l/    1  b^ 


'-a: 


2   r  ar^  "  2  ^9 


(2-10) 
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where  a      ,  o  ^,  a^„,  represent  respectively  the  radial^  shear  and  circum- 
rr   rQ        0&' 

ferential  (hoop)  stresses  In  the  medium. 

In  the  following  development  of  the  equations,  the  incoming  wave, 
considered  a  step  function  of  time,  is  either  a  dilatational  or  shear  wave. 
The  step  wave  was  chosen  as  a  convenient  form  from  which,  through  the  principle 
of  superposition  using  Duhamel's  integral,  the  effects  of  any  time  dependent 
wave  may  be  approximated.   It  is  shown  hy  Kolsl^  (3)  that  for  either  the 
incoming  dilatational  or  shear  wave,  the  reflected  wave  from  a  plane  body  must 
have  both  dilatational  and  shear  components  in  order  to  satisfy  boundary  con- 
ditions.  The  same  holds  in  the  present  case, 

2 . 11  Dilatational  Wave 

The  orientation  of  the  incoming  dilatational  ( compress ional)  wave 
with  respect  to  the  spatial  coordinates  and  time  are  shown  in  Fig.  1.   Owing 
to  symmetry  about  the  x  axis,  the  radial  stress,  hoop  stress,  and  the  radial 
component  of  displacement  are  even  functions  of  9,  while  the  shear  stress  and 
tangential  component  of  displacement  are  odd  functions  of  9. 

In  the  expressions  previously  given,  the  dilatational  potential  may 
now  be  separated  into  two  parts 

Cp  =  cp.   +  cp   , 

^   ^m   ^out 

where 

cp.   =  -  CJ  X  (x+c,t) 

represents  the  potential  of  the  incoming  step  wave,  and  cp    represents  the 
potential  of  the  outgoing  dilatational  wave. 

The  displacement  and  velocity  of  a  particle  behind  the  plane  wave 
front  can  be  derived  from  the  potential  of  the  incoming  wave,  with  notation 
defined  in  Fig.  1 
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u  =  -  a    "/-    ^xi-c,  t) 
X      p     ■   1 

(a-ii) 

P 

X     pc 


The  velocities  in  polar  coordinates  become 


u   =   ■ 

a 

-  -2—  cos 

0 

— *—  sin 

V    = 

(2-12) 


The  radial  and  tangential  displacements  at  the  boundary  r  =  R  due  to  the 
incoming  plane  stress  wave  are  expanded  in  Fourier  series  as  a  function  of 
time  as 


i(R,t,e)  =  y   a^(R,t)  cos  n0 


n=0 

(2-13) 


(R,t,0)  =  y   b^''P.,t)  sin  nQ 


V .  ,  .  . 

n  =  l 


and  the  velocities  as 


L(R,t,e)  =  y   a^(R,t)  cos  n0 


n=0 

(2-1^) 


v(R.t,e)  =     )      i    (R,t)  sin  nQ 


n=l 
The  Fourier  coefficients  are  determined  as  a  function  of  9-,  ;,  a  measure  of  the 
degree  of  envelopment  of  the  shell  by  the  incoming  stress  wave  (Fig.  1),  by 
evaluation  of  the  integrals 


11^ 


2a 


a 


n     npc 


2a 


cos  9   cos  n9  d9 


1  ^0 


1 


— t^     /   sin  9  sin  nQ  d.Q 


n     Jipc^  oq 

The  indicated  integrations  are  performed  using  the  orthogonality  properties 
of  the  sine  and  cosine  functions  to  get  the  following  result 


n(in) 


:tpc- 


n(inj 


rtpc- 


sin  9 


1 

sin  29 


9  +  

1      2 


1 


sin(n-l)9    sin(n+l)0 


n-1 


n+1 


sin  29 


1 


1      2 
sin(n-l)9    sin(n+l)9 


n-1 


n+1 


n  =:  0 


n  =  1 


n  =  2,  3, 


(2-15) 


n-1 


n  =  2,  3, 


The  subscript  (in)  has  been  added  to  designate  the  coefficients  due  to  the 
incident  wave.   9   is  a  time  dependent  variable  which  varies  from  0  to  it 
during  the  transit  of  the  incoming  wave  across  the  cavity^  after  which  it 
becomes  a  constant  equal  to  Jt. 

The  stresses  in  the  medium  behind  the  wave  front  are  as  shown  in 
Figo  1;  where  v ,      ^  ratio  relating  the  stress  in  a  direction  perpendicular  to 
the  direction  of  wave  propagation,  is  derived  from  the  assumption  that  there 
is  no  strain  behind  the  wave  front  in  the  direction  parallel  to  the  wave  front, 
It  is  related  to  Poisson's  ratio,  v^  of  the  medium  in  the  following  way 


*     The  bar  over  the  symbol  does  not  indicate  a  vector  quantity, 


^  "  1  -  V 

The  stresses  behind  the  wave  front  are  in  polar  coordinates 

a        =   -   a   (cos  9  +  V  sin  9) 
rr      p 

cTan  =  -o  (sin  9  +  V  cos  S) 
yy    p 
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(2-16) 


(2-17) 


stresses  in  the  medium  as  a  function  of  radius  and  time  are  now 
expressed  in  terms  of  Fourier  series 


o^^{r,t,e)    =  2^  A^(r,t)  cos  nQ 


n=0 


uu 

a^e(r,t,e)  =    2^  B^(r,t)  sin  n9 


(2-18) 


n=l 


OQQ{r,t,d)    =  ^  C^(r,t)  cos  nQ 


n=0 
and  the  coefficients  determined  in  the  same  manner  as  used  previously  for 
velocities  are 

n  =  0 


(l+v)0^  +  (~^)  sin  29^ 


l-V 


.  V  in )     2Jt 


^  •  <    (l-v)e^  +  (1+v)  sin  29^  +    (-^)  sin  k9 ^ 


n  =  2 


(l-v)sin(n-2)0    (l-v)sin(n+2)0    2(l+v)  sin  n9 
n-2 


(2-19) 


n+2 


V 


n  =  1,5A,5, 


B 


o   (l-v) 
P 


sin  i+S 


1 


k 


n  =  2 


(in)     23t 


sin(n-2)e    sin(n+2)e. 


n-2 


n+2 


n  =  1,3A,5, 
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'n  (  in ) 


a 

_V 

2n 


■  l-v^ 


{l+v)e^   -    (^)  sin  20^ 

^(v-l)e^  +  (l+v)  sin  29^   +  (^)  sin  ke_^ 
(v-l)sin(n-2)0    (v-l)sin(n+2)0^   2(l+v)sin 


n  =  0 


n  =  2 


ne 


1 


n-2 


n+2 


where  tne  subscript  (in)  again  refers  to  the  incident  stress  wave. 

As  noted  earlier^  the  reflected  and  diffracted  wave  is  represented 
as  a  diverging  or  outgoing  cylindrical  wave  from  a  line  source  whose  origin  is 
the  axis  of  the  cylinder,  and  includes  both  dilatational  and  shear  components. 
The  outgoing  potential  functions  are  expressed  as  infinite  series,  each  term 
of  which  will  henceforth  be  called  a  mode 


oo 

cp  _^  =  )   f  (r.t)  cos  n© 
^out    Z_j   n   ^  ' 


n=0 


ijf   ,  =  7   g  (j^jt)  sin  ne 
out    /  ,  °n '    ' 

n=I 


It  is  shown  in  the  appendix  that  the  coefficients  are  of  the  form 

00 

(zl)° 


f  = 

n     n    J 
^1     ° 


F^( 


t  -  ' — •  cosh  u^  )  cosh  nu,  du-, 
c        1'        11 


(-1)^   r   n     r 
^n  "   n    J   '^  ^^  "  ~  ^°^^  ^2-^  ^°^^  ™2  ^^'c 


Cg   ^0 


where  F  (t  -  • —  cosh  u  )  and  G  (t 


(2-20) 


(2-21) 


'1 


1' 


' —  cosh  u  )  represent  the  nth  order 
Cg        2 


derivative  of  the  respective  functions,  n  being  equal  to  the  mode  number 
considered. 


IT 


To  express  the  velocity  components  of  the  medium  at  the  ■boundary- 
due  to  the  outgoing  waves  as  in  equations  (2-1^^-),  one  differentiates 
equations  (2-9)  with  respect  to  time  and  substitutes  the  expressions  given 
above  for  the  potential  functions.   This  results  in  the  following  for  the 
velocity  coefficients 

a  /    X  =  -  ^"  :;    /  F   (I-,)  cosh  u-,  cosh  nu^  du^ 
n(out)      n+1  J  ^^r  1  11 

^1 

00 

+  ■■  "  ■   /   G   dp)  sinh  u  sinh  nu  du 


"2 


(2-22) 


b    /      ^\    =   -    ^"   1.         /     f'^      (|,  )    sinh  u^    sinh  nu^    du 
n(out)  n+1      J_  '1  1  11 

c  0 

00 

^    /    ^'"''^(^2^    "°"^  ^2   "^^'^  ^^2   ^^2 


+ 
c 


2 

where  (^-i)  ^i^d.  (^_)  represent  the  argijments  of  the  respective  functions. 
The  subscript  (out)  refers  to  the  outgoing  wave. 

The  stresses  in  the  medium  due  to  the  outgoing  waves  are  also 
expressed  in  series  form  as  in  equations  (2-l8)^  the  coefficients  being 
found  by  substitution  into  equations  (2-10) 

(-1)  1^ 


n  ( out )    n+2    ^^ 
c,      0 


/      Y)^0  \  p 

/   F   (I-,  )  cosh  nu-,[—  +  2  cosh  u-,  1  du^ 
Jo       1        1m-  l'   1 


-  ^  "njg^  J      ^^'^    ^^2^  ^^"^^  ^2  ^^"^   "^2  ^^2 

(2-23) 


"2 
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B    /         \    =    ^ "    ^ J^      /      f"^      (i^)    sinh  2u,    sinh  nu,    du^ 
n(out)  n+2      J^  ^   1  1  11 

^  c  0 

00 


=2 


oo 

C  ,    .  =  -^^i^'  [    F^'+^C^J  cosh  nn   [-  -  2  sinhV]  du 
n(out)     n+2   J       ^1'        1  ^IJ.  l'   1 

^1 

00 

+  ~^^2      J      ^^'^    ^'^2^    ^^"^^  ^^2  ^^'^^  ^^2  ^^2 
Cg      0 

2.12  Shear  Wave 

The  orientation  of  the  incoming  plane  shear  wave  with  respect  to 
the  spatial  coordinates  and  time  are  shown  in  Fig.  2.   In  this  case  the  shear 
stress  and  tangential  component  of  displacement  must  "be  even  functions  of  9 
while  the  radial  stress,  hoop  stress,  and  radial  component  of  displacement 
are  odd  functions  of  d. 

The  incoming  wave  expressed  in  terms  of  a  shear  potential  function 

is 

A|f.   =  a  X  (x  +  c  t)  (2-21+) 

in    s   s      2  ^  ' 

from  which  the  velocity  of  a  particle  behind  the  wave  front  in  polar 

coordinates  may  be  derived  as 

a 

s 
u  =  — — •  sm  Q 

pCg 

(2-25) 

S         n 
V  =  COS  9 

pcg 

The  stresses  in  the  mediimi  behind  the  wave  front  are  as  shown  in  Figo  2.      In 
terms  of  polar  coordinates  they  are 
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a        =  a     sin  29 
rr    s 


a   ^   =  a     cos  29 
r9    s 

o^„   -   -a     sin  29 
99  s 


(2-26) 


Velocities  and  stresses  due  to  the  incoming  plane  stress  wave  as  functions 
of  radial  distance  and  time  are  expanded  in  terms  of  Fourier  series  in  a 
manner  similar  to  that  of  the  dilatational  wave,  except  that  because  of  the 
difference  in  symmetry _,  the  sine  and  cosine  terms  are  interchanged.   The 
coefficients  are  also  found  in  a  similar  manner  and  are  given  here.   The 


coefficients  for  the  velocity  series  are 

sin  29. 

'i-  — 


X  J 


1 


n(in)   rtpcg   ^ 


sin(n-l)9    sin(n+l)9 


n-1 


n+1 


sin  9 


1 

sin  29 


n(in)   jTpc, 


X  <  9.  + 


1 


sin(n-l)9    sin(n+l)9 
-  + 


n-1 


n+1 


n  =  1 


n  =  2,3, 


n  =  0 


n  =  1 


n  =  2,3,.. 


(2-27) 


and  the  coefficients  for  the  stress  series  are 

sin  k9. 


A 


9 


1 


n(in)    Ti 


1      \ 
sin(n-2)9    sin(n+2)9 


n-2 


n+2 


n  =  2 


n  =  1,3,1^,5,... 


B 


n(in) 


sin  29 


1 


sin  hQ 


1 


1      h 
sin(n-2)9    sin(n+2)9 


n-2 


n-l-2 


n  =  0 


n  =  2 


(2-28) 


n  =  1,3A,5, 
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C  /.  N  =  -  A  /  .  S 
n(in}      n(in; 

9      is  a  time  dependent  variable  similar  to  that  for  the  dilatational  wave 
except  that  now  it  is  a  function  of  the  position  of  the  incoming  shear  wave, 
which  travels  with  velocity  Cp. 

The  outgoing  potential  functions  are  expressed  modewise  as 


CD  ^  =  )   f  (r,t)  sin  nS 
^out    /_,   n   ^ 


n=l 


(2-29) 


^out 


n=0 
where  the  coefficients  are  of  the  form 


00 

f  =  -^ —     /  F  (t  -  ■ —  cosh  U-,  )  cosh  nu^  du^ 

n     n    J    ^    C-,       1         11 


n 

n 

=1 

^0 

OQ 

=1 

'n 

(-1)" 

n 
^2 

r 
"2 

(2-30) 


cosh  u  )  cosh  nu  du 


The  terms  have  the  same  significance  as  in  the  case  of  the  dilatational  wave. 
Modal  coefficients  are  derived  as  before.   The  coefficients  for  the  velocity 
terms  are 


a  /   ,  V  =  -  -^ — ~     /   F   (i-,)    cosh  U-,  cosh  nu,  du, 
n(out)       n+1  J_     ^1^       1        1   ] 
c      0 


00 


(2-31) 


21 


00 

b  /    N  =  ^~'^K       /  'P^^    (^-i)  sinh  u,  sinh  nu,  du 
n(out}     n+1  J  ^   1  1        -I-   J 


c     '^O 
1 


uu 

^"  \      /  g'^'^  (^„)  cosh  u_  cosh  nu_  du^ 
n+1  J„      2        2        ^   £ 


Cg-  -0 


and  the  coefficients  for  the  stress  terms  are 


00 


=1      ° 


CO 

+  ^^^—^     /  G'^'^^d^)  sinh  2u^  sinh  nu^  du, 
n+2   J_       2         2        2   ^ 


=2      ° 


(2-32) 


oo 

/   xn    p 

B  /    \  =  -  '  "  ^ J^      /   F   (i^)    sinh  2u,  sinh  nu-,  du, 
n(out)       n+2   J       ^1'        1        11 

00 

-  \l2      /  ^^'^    ^'^2^  ^°^^  '^2   '^°^^  ^^^2  ^^2 
Cg     0 

00 

C  ,   ^x  =  ^"^)j-   /  f'^"^^(^J  cosh  nuj-  -  2  sinh^  u,  ]  du, 
n(out)     n+2   J„      ^1'        l^\x  l""    ] 

^1 

00 

-  "^ jg^  J   g'^"^  ii^)    sinh  2U2  sinh  nu^  du^ 


2.2  Equations  for  the  Shell 

In  the  following  derivations  of  the  equations  of  motion  for  the 
shell,  it  is  assiimed  that  the  thickness  of  the  shell  is  small  in  relation  to 
its  radius.   This  assumption  permits  description  of  the  behavior  of  the  shell 
in  terms  of  its  middle  surface,  and  effectively  concentrates  all  of  the  mass 
in  a  line  thickness.   This  should  be  kept  in  mind  when  interpreting  the  results 
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The  effects  of  shear  and  rotatory  inertia  are  not  considered,  since  they 
affect  only  the  very  high  modes,  which  are  not  considered  in  this  study. 

The  deformation  of  any  point  (S)  on  the  shell  can  then  be  described 
by  specifying  the  components  of  displacement  (Fig.  3) • 

u  =  radial  displacement,  positive  outwards 

V  =  tangential  displacement,  positive  in  the  direction 
of  increasing  9. 

The  strain  energy  of  the  shell  per  unit  length  can  be  expressed  as 

—     ?jt         —     2it 
E  I   r    o      E  A       n 

u  =  -|-  /  (k)^Rde  +  -|-  /  {e  rme  (2-33) 

where     Es  =  "plane  strain"  modulus  of  elasticity  of  the  shell 

I  =  moment  of  inertia  per  unit  length  of  shell 

A  =  area  of  the  cross  section  per  unit  length  of  shell 

k  :=  change  in  curvature  of  the  cross  section 

e„  =  extensional  strain  in  the  9   direction 

R  =  radius  of  the  cylinder 

Flugge  (2)  gives  the  expressions  for  the  curvature  and  extensional  strain 


terms  as 


-,   d  u 
k  =  -2  (— ^  +  u  ) 


and 


1  /   S       \ 

^9=     R  ^^9-  +  %) 

correct  to  first  order  terms.   The  strain  energy  equation  in  terms  of  the 
components  of  displacement  then  becomes 

iT  T        2Jt     .2  -  ,        2Jt   . 

2R^    ^0        ae  ^  2R      Jq      dy  s 
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The  kinetic  energy  of  the  shell  can  "be  expressed  as 

2Jt 

T  =  I  f      (u^  +  v^)  Rde  (2-35) 

2  Jq   '  s    s 

where 

m  =  mass  of  shell  per  unit  surface  area 

2.21  Dilatational  Wave 

In  view  of  the  symmetries  of  the  problem,  the  displacement  com- 
ponents of  the  shell  can  be  expanded  in  terms  of  the  following  Fourier  series 


u  (9,t)    =  )   a   (t)  cos  nd 


n=0 

(2-36) 


V 

s 


(0,t)  =  y^    b^^(t)  sin  n9 


n=l 


Equations  for  the  strain  and  kinetic  energies  of  the  shell  are  then  expressed 
as  the  following  quadratic  forms  in  the  generalized  displacements  and 
velocities 

— .  _  —  C)0  —  cc 

E  I:r       E  Art        ^         E  In    v-i      ^  ^  ^       E  An 


mm        iii  Art  £,  irt     ^  Ji;  Art     r^ 

-o3  R    '      so  ^^3    Z_j      sn^  2R      ^  sn  sn"^ 

(2-37) 


R''  2R'   ,  , 

n=l  n=l 


^=2    mrtR  \   /o2   ,  ^2  V 
T  =  mrtRa   +  -^r-  /   (a   +  b   ) 
so    2   /_,    sn    sn 

n=l 

The  Lagrangian  equations  of  motion  in  terms  of  the  generalized  displacement 

coefficients  a   and  b   are 
sn      sn 

dt  ^^a   ^  "^  ^a   -  '^^n 
sn      sn 

d  ,aT  .    au   -  (2-38) 

dt  ^ab    ^  ^  w~  -  % 

sn      sn 


2h 


wi-'ere  Q  and  Q  are  generalized  forces  corresponding  to  the  displacement 

n      n 

terms  a   and  b   ,  respectively.   The  following  equations  of  motion  result 
sn      sn' 

E  I   E  A  Q 

/  s      s  \  o                             ^    p. 

a   +  (■— r  +  • — ^)  a   =  - — -                                                           n  =  0 

so   ^  ^k         J2.'  so   2;nnR 
mP,    mR 


(2-39) 


El       2       E^A  Q^ 

a   +  -V  (l-n  )   a   +  -^  (a   +  nb   )  =  —^        n  =  1,  2,  »  .  . 
sn    J^  sn    ^2  ^  sn     sn'   mnR 

mR  mR 

nE  A  Q 

b   +  -~  (a   +  nb   )  =  —^  n  =  1,  2,  .  o  . 

sn    ^2    sn     sn    mitR 
mR 

Now  to  determine  the  generalized  forces  Q  and  Q  ,  the  principle 

n      n 

of  virtual  work  is  applied.   Consider  a  virtual  displacement  corresponding 

to  an  increment  5a   of  coordinate  a 

sn  sn 

Virtual  work,  by  definition  =  Q   °  5a 

'     ^  n     sn 

The  external  forces  acting  on  the  shell  are 


a   =  7  A  cos  n0 
rr    /  ,   n 

n=0 


Q   r^   -     )      B  sin  n9 
xd         /   ,   n 

n=l 

where  A  and  B  represent  the  sum  of  the  forces  on  the  shell  due  to  the 
n      n 

incoming  and  outgoing  stress  waves 

A  =  A  (in)  +  A  (out) 
n    n^   '    j^\    / 

B  =  B  (in)  +  B  (out) 
n    n^   '    ]^\    / 

The  coefficients  A  and  B  have  been  discussed  previously.   The  work  due  to 
n      n  !:■  J 

a  virtual  change  of  the  generalized  coordinate  a   is 

sn 

2n 

Virtual  Work  =   /   A  cos  n9  •  5a    °  RdS 
J  Q   n  sn 
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The  generalized  force  term  is  determined  by  equating  the  two  expressions  of 
work,  from  which  we  get 

Q  =  2nR  A  n  =  0 

o        o 

Q^  =     jrR  A^  n  =  1,    2,    ... 

and  by  a   similar  computation 


Q     =  jtR  B 


n 


The  equations  of  motion  now  become 


E  I     „  2  E  A               A 

a   +  — r  (1-n  )  a   +  — --  (a   +  nb   }  =  — 

sn    T.^  sn    _,2   sn     sn'    m 

mR  mR 


nE  A  B 

b   +  — —   (a   +  nb   }  =  • — 

sn    ^2    sn     sn'  m 
mR 


(2-i^O) 


2.22  Shear  Wave 

Because  of  the  difference  in  symmetry  associated  with  the  shear 
wave  as  compared  to  the  dilatational  wave,  displacement  components  in  the 
case  of  the  incoming  shear  wave  are  written  as 


u 
s 


(0,t)  =  )   a   (t)  sin  nS 
'      l__i      sn 


n=l 

(2-^1) 


Vg(9,t)  =  ^  b^^(t)  cos  n9 


n=0 

Expressions  for  the  strain  and  kinetic  energies  in  terms  of  the  generalized 

coordinates  a   and  b   now  become 
sn      sn 

—  CO  00 

n=l  n=l 

(2-J+2) 
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(2-^5) 


n=l 
The  equations  of  motion  determined  in  a  manner  similar  to  that  of  the 
dilatational  wave  are  then 

E  I     ^  2       E  A  A 

s   /^   2v         s   /  ,   V     n 

a   +  — r  (1-n  )   a   +  — ^  (a  -  nb   )  =  — 

sn    ^k   ^           ^    sn    „2  ^  sn  sn'    m 
mR              mR 

nE  A  B 

Id ^^   (a   -  nb   )  =  — 

sn    „2    sn     sn'    m 
mE 


2o25  Effect  of  Additional  Mass 

The  equations  of  motion  derived  thus  far  assumed  no  mass  within  the 
shell.   The  simplest  way  to  consider  the  effect  of  additional  mass  is  to 
assume  that  it  is  distributed  within  the  shell  symmetrically  about  the  axis. 
The  total  quantity  can  be  assumed  equal  to  23iRm'  ,  with  m'  the  magnitude  of 
the  added  mass  measured  in  terms  of  the  surface  area  of  the  shell.   If  it  is 
assumed  to  move  with  the  same  velocity  as  the  mass  center  of  the  shell,  its 
kinetic  energy  can  be  expressed  as 

2ji  2 


m.  =  2"Rni'   1 

2,2 
k-si 


/   (-U  cos  9  +  V  sin  9)d9 
Jn     s  s 


(2-^1+) 


for  the  case  of  the  incoming  dilatational  wave.   If  u  and  v  are  replaced  by 
their  Fourier  expansions,  and  if  the  indicated  integration  is  carried  out, 
the  kinetic  energy  T'  due  to  the  additional  mass  may  be  represented  as 

4   ^  sl     si  si    sl'^ 

This  kinetic  energy  must  be  added  to  the  kinetic  energy  of  the  shell  which 
was  derived  earlier.   Since  the  above  expression  affects  only  the  n  =  1  mode. 
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the  other  modes  need  not  be  considered ^   With  this  additional  kinetic  energy, 
the  equations  of  motion  for  the  n  =  1  mode  become 

E  A  A, 


E  A  B^ 

mE 


(2-il5) 


If  now  we  add  the  above  equations,  we  get 


2E  A  A^   B^ 

mR 


or 


2E  A  A    B 

raR 


and 


2E  A  A,    B^ 

S      /  T       N          1         1 

b  ,  =  -  a  ^  -  • — -r-  (a  +  b  ^  )  +  -—  +  — 

si      si     „2   si  si'    m    m 
ran. 

By  making  the  appropriate  substitutions  for  a   and  b   ,  the  equations  of 

motion  can  be  written  as 

E  A              A,  ^   ,    B, 

s   /      -,   \    1  2m+m      1     m 

a  -,  +  — —   (a  ^  +  b  -,  )  =  —  •  -TT— — TT  +  ~~  '  7T7~~ — TT 

si    ^2   si    si     m  2(m+m  )    m    2(m+m  ) 
mR 


(2-k6) 


•«     ^s^  .      ^   V   "^1      m'     ^1    2m+m' 


si    „2   si    si  m    2(m+m' )    m    2(m+m' ) 

mR  .    /         \    / 

Note  that  these  equations  are  similar  in  form  to  the  equations  derived  earlier 

without  the  additional  mass;  the  effect  of  the  additional  mass  merely  alters 
the  right-hand  side. 

In  the  case  of  the  incoming  shear  wave,  the  kinetic  energy  of  the 
additional  mass  is 


_  2nRm'  1 

2   ^«2 
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2«  2 

f   (u  sin  9  +  V  cos  9)  dd  (2-i^7) 

0  ^  3  J 


and  by  a  similar  method,  the  equations  of  motion  for  the  n  =  1  mode  are 
written 

V  /  V   ^  _  5.    an+m'    \          ^  m' 

^sl  "^  „2  "^^sl  "   sl^  ~  m    2(m+m'')  "  m    2(m+m' ) 
mE 

(2-1^8) 

..     ^s"^  ,  y,      \    -       tl           m'      \          2m+m' 

"^sl  ■  „2  ^^sl  ■  ^sl''  "  '  m  *  2(m+m')  "^  m  "  2(m+m'  ) 
mR 

The  significance  of  the  additional  mass  on  the  numerical  results  is  discussed 

in  a  later  chapter.  Inclusion  of  a  flexible  support  for  the  additional  mass 

is  also  possible  by  changing  the  above  equations;  however,  this  problem  is 

not  treated  here. 


2.3  Boundary  Conditions 

It  is  assumed  that  the  shell  is  attached  to  the  medium  permitting 
no  differential  displacements  at  the  boundary  between  the  two.   Thus,  con- 
tinuity of  stresses  and  displacements  are  maintained  at  the  boundary.   A 
convenient  way  of  satisfying  the  boundary  condition  of  equal  displacements 
of  the  shell  and  medium  is  to  equate  the  corresponding  velocity  components 
of  the  shell  and  medium.   This  is  done  modewise  for  both  the  dilatational 
and  shear  incoming  waves  by  merely  equating  the  coefficients  of  the  velocity 
terms  as  follows 

a   =a(in)+a  (out) 
sn    n^   '    n^    ' 

(2-^9) 

b   =  b  (in)  +  b  (out) 
sn    n   '    n    '^ 
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2 . h     Summary  of  Equations  in  Non-Dlmensionalized  Form 

It  is  convenient  to  express  the  equations  derived  earlier  in 
non-dimensionalized  form  by  introducing  certain  dimensionless  variables  and 
parameters .   Let 

a     E  tc 

a  =  -1^  •  -^  T  =  ^  (2-50) 

n    R    cf  -K 

P 

b     E 
ft sn  .  _s 

P„  -  R    0^ 

be  the  non-dimensionalized  form  of  the  displacement  and  time  variables. 
Parameters  are  expressed  as 

E 


^E 

E 

^P 

Ps 

P 

\ 

(l+v)(l-2v) 

ii-v; 

\ 

t 

""  R 

\ 

A 

~  R 

\ 

I 
%5 

(2-51) 


The  last  two  parameters  given  above  become  unnecessary  when  considering  an 
unstiffened  shell  since  for  this  case  A  and  I  are  given  by 

A  =  t      and      ^  "  o^  ^^  (2-52) 

Although  the  equations  were  derived  considering  an  unstiffened  shell  of 
uniform  thickness,  the  response  of  a  shell  with  thin,  closely  spaced  stiffeners 
can  also  be  approximated  by  specifying  A  and  I  separately  from  t.   The 
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equations  of  interest  in  non-dimensionalized  form  are 

d  =  PI  (0j  +  Ql  (F,G) 
n     n^  1'     n^ 

P  =  P2  (ej  +  Q2  (F,G) 
^n     n^  1'    ^  n^  '  ' 

a    =  Nl  (a  ,p  )  +  P5  (9)  +   Q3  (f,g) 

S  =  N2^(a  ,8  )  +  Pi|„(0j  +  Q^^(F,G) 


(2-55) 


where  PI  ,  P2  ,  P3  ,  and  '?\     are  functions  whose  values  are  determined 
n'    n    n'       n 

directly  from  the  position  of  the  incident  wave.   Ql  ,  Q2  ,  Q3  ,  and  Qij-  are 

functions  of  the  outgoing  shear  and  dilatational  waves:  and  Nl  and  N2  are 

°   °  n       n 

functions  of  the  displacements.   These  various  functions  in  non-dimensionalized 
form  are  derived  from  equations  given  earlier.   Given  here  are  those  for  the 
case  of  the  incoming  dilatational  stress  wave. 


PI  (0  )  =  .  .^Ll 

n^  1'      TT 


sin 


1 


1 

sin  29 


1 


sin(n-l)0    sin(n+l)9 
-  + 


i,    n-1 


n+1 


n  =  0 


n  =  1 


n  =  2,3. 


P2  (0J 

n  1 


Ve 

TT 


sin  20 


1 


1      2 
;in(n-l)0    sin(n+l)e. 


n-1 


n+1 


n  =  1 


n  ^   2,3,- 


(2-5^) 


P3,(e,) 


w   .  , 


2^p\ 


,    (1-v)  sin  29 
(l+v)Q-^  +  2 


(l-v)9^  +  (l+v)  sin  20^  + 


(l-v)  sin  1+0. 


1  ■      h 

(l-v)  sin  (n-2)0^   (l-v)  sin  (n+2)0 

— +  _ —   -  ■ 

n-2  n+2 

2 (l+v)  sin  n0 


n  =  0 


n  =  2 


n=l,3,i^,5.. 
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P^,(^) 


sin   kd 


1 


sin(n-2)e^        sin(n+2)e- 


n-2 


n+2 


n  =  2 


n   =  1,3,^,5, 


Nl   (a  ,p  ) 

n^    n'    n 


Ve 


-n.  (a     +  np   )   +  il_(l-n    )      a 


N2  (a  ,p  ) 

n^  n'^n 


Q1^(F,G) 


Ve 


T]    (na     +  n   p   ) 

'An  n 


(-l)      /     F        (^,  )    cosh  u     cosh  nu     du 


c     n+1    p 
+    (-1)      •    ( — )  /      G        (^    )    sinh  u     sinh  nu     du^ 


Q2^(F,G) 


■n        /    .  inJ-O 

=  -(-l)      /      F        (^)    sinh  u     sinh  nu     du 


c      n+1    p      ^_^2 
(-1)'^    •    ( — )  /      G        (Ig)    cosh  Ug   cosh  nu^  du. 


Q3^(F.G) 


(-1)"  ^ 


TT 1; /      F        [i^)    cosh  nu^ 


—  +  2   cosh     u, 
.1^  1. 


du. 


(-1)      n^  c^  n+2     p      ^^2 

2(i+v)t]   Ti^    '    ^^^         Jq   ^        ^^2^    "^^^  2U2   sinh  nu^  du. 


Q.\{'P,G) 


i-lf  n. 


V         /  ■    Tnn+2 


2(i^vTyi;  Jo 


F        (^    )    sinh  2u     sinh  nu     du 


2(l+v)ri  T]^    •    ^^^         j^  G        (^2^    =°^h  ^Ug   cosh  nu^  du^ 
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2 • 5  other  Equations  of  Interest 

Dllatatlonal  Wave 

The  hoop  stress  in  the  shell  is  determined  from  the  relation 

where 

1  /   s      \ 

By  substitution  of  the  Fourier  expansion  of  the  displacement  components  the 
hoop  stress  can  be  expressed  modewise  as  a  fraction  of  the  absolute  amplitude 
of  the  incoming  stress  wave,  as 


np  +  a 


cos  n9  (2-55) 


P    n=0 


The  bending  stress  is  determined  from  the  relation 

E  d  ^^u 
s   /    s      \ 
a  .  =  — p  (■ p  +  u  ) 

which  may  also  be  expressed  modewise  as 

-r^  =  d  ^   (n^-1)  a^   COS  nd  (2-56) 

^       n=0 

where  d  is   the  distance  from  the  neutral  axis   of  the  shell  to   its   extreme 

fiber o 

Stresses  within  the  medium  at  any  radius  r  may  be  determined  by 

summation  of  modal  stresses  as  follows 

00     

a     v-1  A 
rr    \    n       „ 
cos  nU 


P    n=0   P 


(2-57) 
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00     — 

r0    \     n    .    -3 
-; r  sm  n0 


^  n=l 


00     — 

COS  no 


a     z_j   cr  i 

where  A  ,  B  ^  and  C  are  coefficients  previously  described  as  generalized 
n^   n       n  jr        ^  o 

stresses  of  the  incoming  and  outgoing  stress  waves. 


Shear  Wave 

The  corresponding  equations  for  the  incoming  shear  wave  are,  for 
the  shell 


%9 


)   (a  -nS  )  sin  nQ 


n=l 


^Sb       .      \        r.        2 


(2-58) 


d  J      (l-n  )  a     sin  n0 


n=l 
and  for  the  stresses  in  the  medium 


a  ^   A 

rr      \    n    .    ^ 
sin  no 


'  s '     '— ^  I  s 
n=l 


00     — 

rB  \  n 


n=0 


sm  nt 


^  cos  n0  (2-59) 


I  s '     'His 
n=l 


where  all  terms  are  as  defined  before, 
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CHAPTER  III 
METHOD  OF  SOLUTION 

3.1  General 

In  this  chapter  are  presented  the  numerical  techniques  used  in 
solving  the  problem.   Since  the  methods  employed  are  similar  in  theory  for 
both  the  incoming  dilatational  and  shear  waves,  only  the  solution  to  the 
incoming  dilatational  (compressional)  wave  will  be  discussed  in  detail  here. 
The  equations  (2-53)  'to  be  solved  consist  of  two  pairs  of  coupled  integro- 
differential  equations  in  the  generalized  coordinates  of  the  shell  and  the 
potential  functions.   The  integral  terms  in  these  equations  contain  elements 
of  the  outgoing  wave  potentials,  both  dilatational  and  shear,  whose  values 
must  be  determined  through  application  of  the  boundary  equations  at  each 
instant  of  time  considered. 

In  the  numerical  solution,  time  is  taken  to  be  the  independent 
variable,  and  is  expressed  in  terms  of  the  half  transit  time  of  the  incoming 
wave  across  the  cavity,  x  =  tc,/R.   Increments  of  time  are  expressed  as  a 
fraction  of  the  half  transit  time;  if  N  is  defined  as  the  number  of  time 
steps  for  the  wave  to  travel  one  radius,  then  i/W  denotes  the  elaspsed  time 
after  the  ith  step  has  been  taken  (Fig.  k) . 

The  position  angle,  9-,  ,  can  be  defined  for  any  step  in  time  as 

9^   =  arc  cos  (l  -  i/N)  (3-l) 

After  total  envelopment,  i.e.,  for  i/N  >  2,  9     becomes  equal  to  tt. 

3-2  Numerical  Integration  of  the  Potential  Functions 

Lamb  (k)   notes  that  if  a  point  source  f(t)dz  is  located  at  the 
point  (x  =  0,  y  =  0,  z  =  z) ,  its  effect  at  a  distance  r  from  the  origin  in 
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the  xy  plane  can  be  represented  by  the  equation 

r2    2" 

cp  =  '^ f  (t  -  ^-^^)    dz 

which  may  be  represented  graphically  as  In  Fig.  5.   To  obtain  the  effect  of 
a  line  source  of  density  F(t)  on  the  z  axis,  spherically  symmetric  point 
sources  with  the  same  variation  in  time  can  be  assumed  to  be  situated  all 
along  the  axis,  and  their  combined  effect  represented  as  an  integral  over  all 
the  point  sources 


vr  jfz^' 

'1    -s/rS+z^ 


V = i  r  ^(*  -  ^)  ^^^  <>^> 


The  limits  of  minus  infinity  to  plus  infinity  are  shown  here  in  order  to 

represent  the  general  case.   It  will  be  shown  that  for  a  disturbance  beginning 

at  some  definite  time,  the  limits  of  the  integral  may  be  taken  as  finite. 

Now  let 

2    2    2/,     2/  2x 
r  +z   =r(l  +  z/rj 

=  r  cosh  u 
and  the  above  Integral  can  be  written  as 

so 

cp  =     F{t  -  —   cosh  u  )  du  (3-3) 

from  which  the  graphical  Interpretation  in  Fig.  6  is  made. 

Consider  for  the  present  the  determination  of  the  Integral  values 
at  the  boundary  r  =  R.  We  define  t  =  0  to  be  that  time  when  the  first  effects 
of  the  outgoing  stress  waves  (i.e.,  from  the  point  source  at  z  =  O)  reach  the 
boundary.  As  time  increases  the  effects  of  other  point  sources  progressively 
further  away  along  the  z  axis  will  reach  the  boundary  in  the  xy  plane.   Thus 
it  can  be  seen  that  though  the  limits  of  infinity  are  given  in  the  integrals, 
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for  finite  times  the  limits  can  he   represented  by  finite  values.   The  maximiim 

value  of  the  variable  R  cosh  u,  in  equation  (j-^)  that  need  be  considered  can 

therefore  be  represented  as  a  function  of  time.  At  t  ==  0 

R  cosh  u,   =  R 
Im 

that  is,  U-,   is  zero  for  time  zero.   For  any  time  later 
im 

R  cosh  U-,   =  R  +  tc, 
Im         1 

in  the  case  of  the  outgoing  dilatational  wave  (Fig.  6).   R  cosh  u   locates 
the  point  farthest  along  the  z  axis  whose  wave  front  has  just  reached  the 
boundary  in  the  xy  plane,  c,  being  the  velocity  of  wave  propagation.   In 
non-dimensionalized  form  this  can  be  written 

cosh  U-,   =l  +  'r=l  +  — 
Im  N 

where  the  terms  are  as  previously  defined.   Since  the  integration  is  over  all 

point  sources  whose  wave  front  has  reached  or  passed  the  boundary  in  the  xy 

plane,  the  integral  expression  of  equation  (3-3)  can  be  written  with  new 

limits 

^Im 

CD  =   /    F(t cosh  U-,  )  du^ 

^0         "l       ^    ^ 

Let 

cosh  u,  =  1  +  ^  (3-^) 

where  ^  now  becomes  the  variable  of  integration  varying  in  value  from 
0  to  T .   In  terms  of  ^  the  integral  expression  becomes 

,  -   f  pl^  (J.5) 

Note  that  a  singularity  of  the  integrand  occurs  at  ^  =  0,  for  all  values  of  t. 
Paul  (6)  has  presented  a  numerical  technique  for  evaluating  integrals  of  this 


57 


general  type  taking  into  account  the  singular  point.   The  method  is  essentially 
a  modified  trapezoidal  rule  in  which  weighting  coefficients  are  obtained.   If 
the  range  (0,t)  of  the  integral  is  subdivided  into  subranges  each  of  equal 
length  At,  the  above  integral  for  any  time  t  =  IAt  can  be  represented  as 

(i-l)   (m+l)AT 


9 


m^O  mAT 


(3-6) 


where  f(0  is  a  polynomial  of  first  degree  in  ^  which  approximates  the 
quotient 

F(-l+T-0 

n/ 1+^/2 

as  a  linear  function  over  each  subrange.   The  results  of  the  integration  of 
equation  (5-6)  are  then  put  in  the  form 

(1-1) 


9 


AM   •  R   •  F.    +  BM   •  R   ,  •  F. 
m    m    i-m     m    m+1    i-m+1 


(5-7) 


m-O 


where 


AM  =  J^  i2(m+l) 
m    V  2   '  ^    ' 


(m+1) 


1/2    1/2 


m 


At   2 
BM  =  J-—  ^"^ 
m    V  2 


(m+l)5/2-m5/2 


-2m 


/   ,x3/2  5/2 
J    ^  L 

^z   ,n1/2  1/2 
(m+1)  ^  -m  ■ 


(5»8) 


represent  a  set  of  weighting  functions  which  can  be  computed  for  any  interval 
being  considered.   For  large  values  of  m  these  functions  were  written  in  a 
form  more  convenient  for  computational  purposes  by  asymptotic  expansion,  to 
get  the  following  result 

1 


AM 


m 


m 


mAT 
2 

mAT 
2 


_2m 

\_ 
2m 


1 


1 


2      ^  "  Z      \ 
12m    52m'^   6i«n  -I 


2  ^  _  5 


1 


h 


6m    32m    16m  -^ 


(5=9) 
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For  an  illustration  of  the  variation  of  AM  and  BA   functions  with  m,  see 
Fig.  10,  R  is  termed  the  multiplying  factor  which  for  this  particular 
integral  is  equal  to 

All  of  the  various  integral  expressions  for  the  potential  functions  were 
simplified  to  the  form  of  equation  (3-6)  thereby  permitting  use  of  the  above 
weighting  factors  for  all  integrals.   However,  the  expression  for  the  multi- 
plying factor  R  varies  for  the  different  integrals. 

The  integration  for  the  outgoing  shear  potential  ijf  is  handled  in  a 
similar  manner,  except  for  a  minor  modification  which  results  from  the  fact 
that  the  shear  wave  travels  outward  with  velocity  c  . 

Consider  for  illustrative  purposes  some  shear  potential  represented 
as 

00 

il;  =   /   G(t  -  —  cosh  u^)  cosh  u^  du^ 
Jq       Cg       2        2   2 

Since  the  incident  wave  strikes  the  shell  at  time  t  =  0^  the  shear  wave  front 
must  also  reach  the  boundary  at  this  instant  in  ordei-  to  satisfy  boundary- 
conditions.   Thus  at  t  =  0 

R  cosh  u,^  =  R 

2m 

and  for  any  time  later 

R  cosh  u^  =  R  +  tco 
2m         ^ 

using  the  same  reasoning  as  for  the  outgoing  dilatational  wave.   We  can  express 

the  above  equation  in  non-dimensional ized  form  as 

cosh  u_.  =  1  +  k  T 
2m        c 
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where 

is  a  constant  relating  the  velocities. 

Again  since  u  is  a  variable  of  integration  we  can  write 

cosh  u_  =  1  +  k  C 
2        c 

where  ^  as  in  the  case  of  the  dilatational  wave  is  the  new  variable  of 

integration  with  limits  0  to  T .   In  terms  of  this  variable  the  integral  is 

now  written  as 


.^  G(-  i+  X  -  ON/k'e  (1  +  \^^   ^^ 

(3-10) 


^i^  =     1£ 


M     ,/,  ,  !£i 


Note  that  the  shear  wave  starts  from  the  z  axis  at  time  t  =-l/k  which  is 

'  c 

earlier  than  for  the  dilatational  wave  since  it  travels  at  a  slower  speed 
(with  fictitious  material  within  the  shell)  and  must  reach  the  boundary  at 
time  equal  to  zero.   The  numerical  integration  of  equation  (3-IO)  is 
accomplished  in  a  similar  method  to  that  which  was  done  for  the  dilatational 
potential.   The  result  can  be  expressed  as 
(i-1) 


^  = 

m=0 


AM   •  Q   •  G.    +  M   •  Q    •  G. 

m    in    i-m     m    in+1    i-m+1 


(3-11) 


where  AM  and  BM  are  weighting  factors  identical  to  those  derived  previously, 
mm  -t-        ^  7 


and  Q  for  this  particular  integral  is 

-JY    '    (1  +  k  O 


\  = 


cm 


k  C 

1  +  "^ 


2 

The  expression  for  Q  is  again  dependent  on  the  form  of  each  integral  considered. 


1+0 

For  an  illustration  of  the  variation  of  the  F  and  G  functions  with  time, 
see  Figs.  11  and  12. 

3.5  Solution  of  the  Basic  Equations 

The  equations  of  motion  for  the  shell  are  solved  by  an  iterative 
procedure  known  as  the  Newmark  Beta  Method  (5)  with  which  values  for  the 
potential  functions,  and  the  accelerations,  velocities,  and  displacements 
of  the  shell  are  determined. 

In  general  the  method  consists  of  using  a  step-by-step  integration 
technique  over  successive  time  intervals  assuming  a  specific  variation  of 
acceleration  during  each  interval.   If  we  assume  that  at  time  t  =  (i-l)Ax  all 
values  of  potential  functions,  accelerations,  velocities  and  displacements 
are  known,  the  method  becomes  that  of  determining  the  corresponding  values 
at  a  time  t  =  IAt .   For  an  assumed  linear  variation  of  acceleration  over  each 
interval  the  equations  of  interest  are 

At 


.   =  a     .  ,  +  -—  (a  .  ,  +  a  . ) 
n,i    n,i-l    2    n,i-l    n,i 

.  =  a     .  -,  +  Ata  .  ^  +  -7=-  (2a  .  ^  +  a  . ) 
n,i    n,i-l      n,i-l    5     n^i-1    n,i 


•       •        At 

p   .=p   .^+—  (p   ..+P   .) 

n,±       ^n,i-l    2    n,i-l    n,±' 


(3-12) 


At^   .. 

P.  .  =  P.  .„n  +  ^-^P.  ._n  +    -T-    (2P.  .  n  +  P.  ^) 


3   .  ^  +  AtP   .  ,  +  —r-   (2p   .  ,  +  p 

n,i-l     "^n^i-l    6     n,i-l    n. 


where  n  denotes  the  mode  and  i  the  time.   Note  that  in  the  above  equations 

the  velocities  and  displacements  are  expressed  in  terms  of  a   .  ^  and  B   .  ^  , 

n,i-l     *^n,:L-l' 

acceleration  components  which  generally  are  unknown. 

Values  for  the  acceleration  components  at  time  t  =  iAt  are  assumed 
and  the  velocities  and  displacements  determined  from  the  above  equations;  then 
by  use  of  the  boundary  equations  to  determine  the  values  of  the  potential 


1^1 

fiinctionS;  values  of  the  acceleration  components  can  te  computed  from  the 
equations  of  motion.   The  computed  values  are  then  compared  with  the  assumed 
values  and  if  they  agree  or  are  within  an  arbitrarily  established  tolerance, 
the  step  is  completed^  if  not,  the  values  of  'd      .    and  3   .  just  computed  are 
taken  as  new  assumed  values  and  the  cycle  repeated  until  the  criterion  is 
satisfied.  A  detailed  discussion  of  the  procedure  follows. 

The  equations  of  continuity  at  the  boundary  at  time  t  =  IAt  are 
given  in  equations  (2-4-9)  ^-i^d  are  rewritten  here  in  the  form 


.  =  PI  (9J  +   /   Ql  (F_,F^,...F.;  G„,G,  ,  ...G.)  d^ 
n,i     n^  1'    J  n^  0   1^     i   0  1      i 


0 

(>15) 
iAT 

P   .  =  P2  (e,  )  +   /   0,2   (F_,F,  ,...F.;  G„,G,  ,...G.)  d^ 

0 

where  PI  and  P2  denote  functions  related  to  the  particle  velocities  in  the 
n       n 

medium  at  the  boundary  due  to  the  incoming  stress  wave;  their  values  may  be 
determined  directly  at  any  time.   The  integral  expressions  denote  the  con- 
tribution due  to  the  outgoing  waves .  All  values  of  F  and  G  are  known  except 
F.  and  G.  which  are  to  be  determined  from  the  above  equations.   These  equations 
are  now  written  isolating  the  unknown  values  of  F  and  G  terms 

iAT 

d   .  =  PI  (ej  +   /   Ql  (F_,F^,...F.  ■,    G„,G^,...G.  J  d^ 
n,i     n^  1     J    ^  n  0^  1'     i-l   0^  1'  i-l 


At 
At 


(3-U) 


/  Ql  (F.  ^,F.;  G.  ^,G.)  d^ 
J  (^        n   1-1   1    1-1   1 


0 

IAt 

P   .  =  P2  (ej  +   /   Q2  (F_,F^,,.oF.  ,;  G^,G,  ,...G.  J  d^ 
•^n^i     n^  1'    J  n  O'^  1^     i-l   0^  1^    i-l' 


^At 
At 

+ 


/   Q2  (F.  ^,F.;  G.  ,,G.)  d^ 
Jq    n^  1-1^  x'      1-1^  1- 
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The  Integrals  with  limits  At  to  iAx  are  evaluated  using  the  numerical 

technique  described  in  a  previous  section. 

The  next  step  is  to  assume  values  for  <DL      .    and  B   .  ,  the  radial  and 

n^ 1      n, I 

tangential  components  of  acceleration  of  the  shell,  from  which  velocity  and 

displacement  components  are  determined  using  equations  (3-I2) .   With  these 

values  known,  equations  (3-1^)  become  effectively  two  equations  with  the  two 

unknowns,  F.  and  G. ,  which  can  then  be  determined. 
'      1      1 

The  equations  of  motion  are  given  in  equations  (2-53)  ^i^d  are  here 
written  in  the  form 

iAx 

°  (3-15) 

iAx 

p   .  =  N2(a   .;  p  .)   +  Pk   (9)   +      /   Qif  (F^,F.  ,...F.;  G^,G.  ,...G.)  d^ 
^n,i       n,i  ^n,x  n  ■   1  J„    n  o^  1'    1'   o'  1'    1'   ^ 

from  which  the  acceleration  components  a      .   and  6   .  are  now  computed.   If 

n,i     ^n,i 

the  computed  values  agree  or  are  within  an  arbitrarily  established  limit  of 

the  assumed  values,  the  step  is  completed;  otherwise  the  computed  values  are 

used  as  the  assumed  values  for  the  next  cycle  of  iteration o 

The  procedure  described  above  requires  that  all  parameters  at  x  =^  0 

be  known,  including  the  initial  values  F  and  G  ■>      Initial  values  of  velocity 
^         ^  00 

and  displacement  of  the  shell  are  specified  to  be  zero.   From  the  short- time 
approximation  presented  in  Section  3-7  initial  values  of  acceleration  and  G 
were  determined  to  also  be  zero,  and  the  values  of  F  to  be  the  following 


F  = 

O      IT 


Ve   J  -1  n  .  0 


(3-16) 


UC-l)""^^  n  =  1,2,... 


The  short-time  approximation  also  indicated  that  the  radial  acceleration 

1/2  3/2 

components  near  x  =  0  varied  as  x  '   and  the  tangential  components  as  x  '  ; 


h^ 


thus  to  improve  the  accuracy  of  the  calculations,  equations  (3-12)  for  the 
first  step  in  time  were  modified  to 


^n,l  =1  ""n.l  K,l   =  5  ^n,l 


2  •  2  • 

n,l   5  n^l  ^n,l   7  "^n^l 


(3-17) 


Without  this  modification,  the  results  for  the  first  step  cannot  be  brought 
into  acceptable  agreem.ent  with  the  short- time  approximation. 

The  n  =  0  mode  is  simplified  somewhat  because  it  has  no  tangential 
component  of  displacement,  and  contains  only  the  dilatational  component  of 
the  potential  functions . 

3.  h     Stresses  of  the  Shell  and  Medium 

Solution  of  the  equations  of  motion  yields  modal  values  of  the 
functions  F  and  G,  and  the  modal  acceleration,  velocity  and  displacement 
components  of  the  shell.   The  stresses  in  the  shell  are  determined  directly 
from  the  displacements  using  equations  (2-55)  ^nd  (2-56). 

Stresses  in  the  medium  at  any  radius  and  time  are  determined  using 
equations  (2-57) •   They  are  here  written  in  the  form 


IV  UU  00 


n=0 


0  ^0 


(>18) 


where  S  (9  )  represents  the  stress  due  to  the  incoming  wave  and  the  integral 
terms  that  due  to  the  outgoing  waves.   All  are  functions  of  both  radius 
and  time.   Consider  the  stresses  in  the  medium  for  any  radius  equal  to  r  at 
time  -;  =  iAx  =  i/N. 

d      (Fig.  7)  is  now  determined  from  the  equation 

Q     -  arc  cos  -  (l  -  i/w)  (3-19) 


kk 


and  takes  on  values  from 

9,  =  -      for  T  =  0 
1   r 

to 

> 


0^  =  TT      for  T  ^  (1  +  -) 


The  integral  terms  must  be  recomputed  for  each  radius  and  time  considered. 
As  an  example  consider  the  integral 


00 

cp  =   /   F(t  -  —  cosh  U-,  )  du. 


and  let  t  =  r—  -  1,  where  k  =  ^/^,    represent  the  nondimensionalized  time 
k  r 

r 

after  t  =  0  required  for  the  outgoing  dilatational  wave  front  to  reach  the 

radius  r  "being  considered.   Then  for  t  <  t  the  integral  is  equal  to  zero, 

and  for  t  >  x   the  maximum  value  offr  cosh  u, )  can  be  written 
r  ^        1 

r  cosh  u^   =  r  +  t '  C-, 

Im         1 


or 

cosh  U-,   =  1  +  k  t' 
Im        r 

where  t'  =  t-t   represents  the  time  after  the  wave  front  has  reached  the 
r   ^ 


radius  r.   We  express  cosh  u-,  as 

cosh  u^  =  1  +  k  C 
1        r 

where  ^  is  the  same  variable  of  integration  as  was  used  previously.   The 

integral  is  then  expressed  in  terms  of  ^  as 

The  numerical  method  of  integration  used  is  identical  to  that  described 


i^5 


earlier  for  the  integral  at  the  boundary^  with  the  result  that  the  same 
weighting  factors  AM  and  EM  are  applicable.   Note,  however^  that  the  limits 
of  the  integral  are  now  0  to  t';  also  the  multiplying  factor  R  for  this 
particular  integral  becomes 

R      = 

m 


r  m 
By  a   similar  analysis^    the    integral  of  the   shear  potential 


ilf  =      /      G(t   -  —  cosh  u^)    cosh  u^  du^ 
Jq  Cg  2  2        2 


can  be  put   in  the   form 


T'      G(      -  j^  +  T-    -   0  VkT    (1  +  k^k^Od^ 


c 


^  =      j        ^ ^ ,  (3-21) 

0  , k  k  C 

'  r  c 


-/^   \|l  + 


rs 


where 

T  '      =    T     -    T 

For   the  potential  ijf 

^s  =r^r-^^ 

c       r 

The  weighting  factors  AM     and  BM     remain  the   same  but   the  multiplying  factor 
DO  mm  j^  J'      to 


0     for  this  particular   integral  becomes 

s/k  k     (l  +  k  k  ^    ) 
Q     _         r   c      r  cm 

/     ^     rem 
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3.5  Time  Dependent  Stress  Wave 

The  problem  considered  thus  far  has  dealt  with  an  incoming  stress 
wave  with  a  step  distribution  in  time.   The  results  obtained  for  the  step 
wave  can  be  used  through  the  application  of  Duhamel's  integral  to  find  values 
for  the  response  of  a  shell  and  medium  to  incoming  waves  with  any  time 
variation.  As  an  example,  a  stress  wave  which  decays  exponentially  with  time 
(Fig.  8)  according  to  the  following  equation  is  considered 

a  (O  =  a   (1  -  f )  e   '°  (5-22) 

•^  o 

T   represents  the  time  at  which  the  stress  wave  decays  to  zero  and  k  is  a 
parameter  which  is  related  to  the  shape  of  the  curve.   Stresses  at  any  time 
equal  to  iAt  can  be  found  by  the  application  of  Duhamel's  integral,  here 
written  in  the  form 

IAt 


oiiA-x)    =  o^^ilAx)   +    J  g^    •  a^^(iAT-T)  dT        (3-23) 


a   is  defined  to  be  that  stress  resulting  from  an  incoming  step  wave  of 

S  b 

amplitude  a        in  equation  (3-22).   In  the  numerical  analysis,  the  time 
dependent  wave  is  approximated  by  a  series  of  rectangular  sections  as  illu- 
strated in  Fig.  8. 


3.6  Description  of  the  Computer  Program 

The  problem  was  programmed  modewise  for  a  high  speed  digital 
computer  (CDC  l60k)   using  Fortran  language,  and  considering  only  the  first 
three  modes.   Input  data  consist  of  the  following  parameters:  -^ 

(1)  The  time  intervals  at  which  computations  are  to  be  performed, 
expressed  as  the  number  of  intervals  required  for  the  incoming  wave  to  travel 


^7 

one  radius  (one-half  transit  time).   The  degree  of  accuracy  achieved  is 
dependent  on  this  parameter,  the  smaller  the  interval  the  more  accurate 
the  results.   However,  the  machine  time  required  for  a  given  number  of  transits 
of  the  incoming  wave  varies  inversely  as  the  square  of  the  interval  size  ap- 
proximately, thus  some  sacrifice  in  accuracy  is  necessary  to  reduce  the  time 
of  computations  required. 

(2)   The  total  time  over  which  the  computations  are  to  "be  performed, 
expressed  as  the  total  number  of  time  intervals  to  be  considered.   Generally 
speaking,  all  values  seemed  to  have  reached  their  asymptotic  (static)  values 
within  ten  transit  times  of  the  incoming  wave  across  the  cavity. 

(5)   The  ratio  of  the  moduli  of  elasticity,  E  /E,  where  E   is  the 

o  S 

"plane  strain"  modulus  for  the  shell  and  E  the  modulus  of  elasticity  of  the 
medium. 

(k)     The  mass  ratio,  p  /p,  where  p   is  the  mass  of  the  shell  and  p 
the  mass  of  the  medium  per  unit  volume. 

(5)  Poisson's  ratio  of  the  medium. 

(6)  The  ratio  of  the  thickness  of  the  shell  to  its  radius.   When 
considering  a  shell  whose  area.  A,  and  moment  of  inertia,  I,  are  not  directly 
related  to  the  thickness,  A  and  I  must  be  specified  separately. 

(7)  The  amount  of  additional  mass  within  the  shell  expressed  as  a 
fraction  of  the  mass  per  unit  surface  area  of  the  shell. 

(8)  The  number  of  radii  to  which  stresses  in  the  medium  are 
desired. 

(9)  The  time  intervals  at  which  output  data  are  desired 0 

(10)  The  angular  increment  at  which  output  data  is  to  be  computed. 
Because  of  symmetry  only  values  between  0  and  I80  degrees  need  be  considered. 


kS 


Additional  input  quantities  in  the  case  of  the  exponentially 
decaying  stress  wave  are  t  and  k^  parameters  which  define  the  duration  and 
shape  of  the  stress  pulse ,  respectively. 

Output  data  consist  of  the  acceleration,  velocity,  displacement, 
and  stress  components  of  the  shell  for  specified  angles  and  times;  and  the 
stresses  in  the  medium  for  specified  radii,  angles  and  times. 


3.7  Short -Time  Approximation 

As  a  check  on  the  accuracy  of  the  machine  solution  for  short  times, 
the  boundary  equations  and  the  equations  of  motion  were  solved  approximately 
by  making  a  series  expansion  of  all  pertinent  functions  in  terms  of  time  as 
the  independent  variable.  Although  the  following  discussion  is  limited  to  the 
case  of  the  incoming  dilatational  wave,  the  basic  principles  are  the  same  for 
either  of  the  types  of  wave  considered. 

The  expressions  for  the  velocities  and  stresses  in  the  medium 
around  the  boundary  due  to  the  incident  wave  all  have  been  written  thus  far 
in  terms  of  9    ,    the  position  angle  of  the  wave.   In  terms  of  nondimensionalized 
time  T,  all  functions  of  d-.    can  be  written  in  terms  of  1   using  the  following 
relations 

cos  0=1-1: 


sin  0     =    V2t 


1      12 


{>2k) 


9^     =    n/2t 


1      32 


The  integral  Aralues  representing  the  effects  of  the  outgoing  waves 
are  also  expressed  as  functions  of  T.   The  following  example  will  illustrate 
the  technique  used  to  accomplish  this.   For  example,  consider  the  transforma- 
tion of  the  integral 


^9 


so 

cp  =   /   F(t  -  —  cosh  u  )  du 


As  was  shown  in  Section  3 '2;  cosh  u   in  the  integral  varies  in  value  for  any 
given  time  T  from  1  to  1  +  t .   It  is  convenient  for  purposes  of  analyzing  the 
integral  to  represent  the  variable  cosh  u  as  follows 

cosh  U-,  =  1  +  UJT  (3-25) 

where  1   is  now  a  fixed  value  in  the  integration  and  u)  is  defined  as  the 
variable  of  integration  whose  value  ranges  from  0  to  1.   The  function 

F(t  -  —  cosh  u  )  =  F(-1  +  t(1-u))) 
^1 

is  expanded  in  terms  of  a  power  series  as 


F(-l  +  t(1-u)))  =  F^  +  2^  7.T  (1-u))- 

i=l 


(5-26) 


where  7.  are  unknown  coefficients  of  the  series.   Since 
1 


du.  = 


n/2t" 


'"  vu  - 


dw 


the  integral  can  now  be  written  as 

1 

cp   =  ^     /       —     F     +  7v  X  (l-w)   +    .  .  . 
\/2    '^0    n/co    -  °  J 


1    --^-H     .. 


doj 


from  which  after  performing  the  indicated  integration 


cp  =  v  2t 


F  +  t(-J  7,  -  ^  F  )  +  t""  ( 
o     3  1   o  o 


(3-27) 


The  integrals  which  represent  the  effects  of  the  outgoing  shear 
wave  can  also  be  transformed  in  somewhat  similar  manner.   Consider  for  example 
the  integral 


50 


^  =   /   G(t 


—  cosh  u  )  cosh  u  du 

ii-\  ^  £_      ^ 


In  this  case  cosh  Up  varies  from  1  to  1  +  k  t  where  k 
Therefore  we  write 


'2   ,     ^'l 
—  and  T  =  ^r— , 
c  R 


cosh  u^  =  1  +  u)k  T 
2         c 


(3-28) 


where  w  is  now  the  variable  of  integration  ranging  in  value  from  0  to  1  in 
the  integral.   The  function 


G(t cosh  u^)  =  G(-  —  +  t(l-oj)) 

2  c 


is  expanded  in  terms  of  a  power  series  as 


G(-  ^+  t(l-u3))  =  G^  +  ^  e.x'd-co)^ 
^  i=l 

where  €.  are  unknown  coefficients  of  the  series.   From  equation  (3-28) 


du. 


N/2k  T 
c 


V  U) 


doa 


The  integral  written  in  terms  of  the  variable  oo  is  now 


^   = 


vkT   pi  , 

c     /     1 


G  +  e,T:(l-uj)  +  , 
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and  performing  the  integration 
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(3-29) 


The  same  basic  technique  is  applied  to  all  the  integrals  so  that  now  the 
effects  of  both  the  incoming  and  outgoing  waves  can  be  written  in  terms  of  x. 
Substitution  into  the  continuity  equations  and  the  equations  of  motion  yield 
the  following  for  any  mode 
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lere  fl.,  f2.,  f3-j  ^nd  f^.  are  coefficients  which  contain  certain  elements 
1^        i'    1        1 

,  G  ,  7.  and  €.  of  the  potential  functions:  il  and  i2  are  known  functions 

)'   o    1        1 

'  the  displacements. 

To  solve  the  above  equations,  the  displacement  components  of  the 
lell  are  expressed  as  Fro"benius  (8)  type  series 
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lere  p.;  q.,  and  c  are  unknown  coefficients.   Substitution  into  equations 
j-30)  results  in  the  following  set  of  equations 
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p.  .  (c+i)-(c+i-l)  .  T^^^-^  +  2,  n(p.,q.)T^-^^  =  ^  f3.(F,G)T'  "  2 
i=l                           1=1               1=1 

00                                                      00                               oo  .      1 

q.  •  (c+l).(c+l-l)  •  t"-'^-2+  ^  i2(p.,q.)T"-'^  =  ^  fi^.(F,G)T^   2 
1=1                             1=1                1=1 


The  coefficient  c  Is  now  determined  by  Inspection.   Then  through  a  step  by 
step  process  which  Involves  the  equating  of  coefficients  of  like  powers  of  i , 
values  of  p.,  q.^  7.,    and  e.  are  determined.   The  following  equations  are  then 
used  to  find  values  of  the  potential  functions,  and  the  displacement  components 
of  the  shell. 


F(t)  ^  F     +     )      7.T^ 
1=1 

00 

G(x)  -%^Y.    '±'' 

1=1 

a(T)  =  2^  Pi  '^'^^^ 

1=1 

(X) 

1=1 

Velocity  and  acceleration  components  may  be  determined  by  differentiation  of 
the  above  expressions  for  the  displacements. 
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CHAPTER  IV 

,■:•    DISCUSSION  OF  RESULTS 

k-.l     General 

Results  of  computations  performed  to  determine  the  effect  of  the 
various  parameters  are  discussed  in  this  chapter. 

Although  equations  presented  throughout  the  study  have  "been  written 
to  include  an  infinite  number  of  modes,  the  greater  part  of  the  actual  calcu- 
lations performed  and  presented  here  are  the  results  obtained  considering  only 
the  modes  n  =  0,  1,  and  2.   It  is  important  to  note  that  during  envelopment  of 
the  shell  by  the  plane  stress  wave,  a  Fourier  series  representation  of  the 
incoming  wave  is  objectionable  in  that  the  series  at  this  stage  is  slowly  con- 
vergent, thus  necessitating  a  large  number  of  modes  to  accurately  represent  the 
plane  wave.   However,  after  passage  of  the  wave  across  the  cavity,  the  Fourier 
expansion  of  the  incoming  stresses  around  the  boundary  results  in  coefficients 
of  all  modes  except  n  =  0  and  2  becoming  identically  equal  to  zero  for  the 
plane  dilatational  wave,  and  coefficients  of  all  modes  except  n  =  2  becoming 
identically  equal  to  zero  for  the  plane  shear  wave.   Therefore,  stresses  due 
to  the  outgoing  waves  in  modes  corresponding  to  those  of  the  incoming  wave 
whose  coefficients  become  zero  must  also  eventually  vanish  at  long  times.   The 
limited  study  conducted  for  modes  greater  than  n  =  2  indicated  that  the 
maximum  effect  of  the  higher  modes  occurs  within  one  transit  time  of  the 
incident  wave  and  rapidly  decays,  thus  contributing  relatively  little  to  the 
maximum  response  of  the  shell  which  occurs  after  several  transit  times. 
However,  for  determining  the  early  time  response  of  the  shell  and  medium, 
the  higher  modes  are  significant  and  should  be  considered  in  further  extension 
of  this  work. 
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In  the  tables  and  figures  to  be  discussed,  quantities  given  in 

non-dimensionalized  units  are  defined  "by  equations  (2-56)  and  (2-5O) .   Stresses 

are  given  in  units  of  the  absolute  value  of  the  amplitude  {\o    |  or  \o    \)    of  the 

P       s 

incident  wave;  a  negative  stress  means  a  compressive  response  to  an  incoming 
compressional  (Fig.  l)  or  a  positive  shear  wave  (Fig.  2).   The  physical 
properties  of  the  shell  relative  to  those  of  the  medium,  as  well  as  the  nature 
of  the  incoming  wave  are  indicated  on  the  graphs.   Unless  otherwise  stated, 
the  shell  is  considered  to  be  an  unstiffened  one  so  that  its  cross  sectional 
area  and  moment  of  inertia  are  related  to  the  thickness  as  given  by  equations 
(2-52).   Except  where  indicated,  there  is  assumed  to  be  no  additional  mass 
within  the  shell.   Numeral  subscripts  denote  the  mode  number. 

The  shell  and  medium  have  been  assumed  to  exhibit  linearly  elastic 
behavior  throughout  their  stress  histories,  which  for  the  practical  problem 
does  not  permit  evaluation  of  any  spalling  or  non-elastic  effects. 

Values  of  stresses  given  are  in  addition  to  those  which  exist  prior 
to  the  arrival  of  the  incident  wave.   For  the  elastic  case,  the  effect  of 
prior  stresses  such  as  those  resulting  from  the  overburden  may  be  taken  into 
account  by  merely  adding  them  to  stresses  caused  by  the  incident  wave. 

For  clarity  of  presentation  and  because  of  the  impracticability  of 
including  solutions  for  all  possible  permutations  of  the  parameters  involved, 
the  discussion  in  this  chapter  is  limited  to  a  few  representative  cases c 

4,2  Modal  Response  of  the  Shell  and  Medium 

Figures  11  and  12  illustrate  the  shape  of  the  modal  components  of 
the  dilatational  and  shear  potentials  obtained  in  the  solution  to  a  typical 
problem.   It  appears  that  a  singularity  occurs  at  one  transit  time  in  the 
case  of  the  F  functions  resulting  in  the  slight  irregularity  of  the  curves  at 
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this  point.   Since  computed  values  of  stresses  were  determined  to  "be  rather 
insensitive  to  relatively  large  variations  in  values  of  the  F  and  G  functions, 
the  effect  of  the  irregularity  would  seem  to  be  slight. 

Figures  IJ  through  20  show  modal  acceleration,  velocity,  displace- 
ment and  stress  components  for  the  shell  and  stress  components  for  the  medium 
at  the  "boundary,  as  they  vary  with  time.   Static  values  shown  were  computed 
using  the  method  given  in  Appendix  B. 

The  high  accelerations  computed  near  the  beginning  are  not  truly 
representative  of  the  actual  case,  since  they  are  the  result  of  assumptions 
made  earlier  in  deriving  the  equations  for  the  shell.   The  shell  was  repre- 
sented by  a  line  describing  its  middle  surface  which  permits  no  variation  in 
accelerations,  velocities,  and  displacements  of  particles  through  the  actual 
thickness.  Also,  no  provision  was  made  for  refraction  of  the  incident  wave 
through  the  shell  lining.   These  limitations  restrict  the  applicability  of 
the  solutions  to  a  shell  whose  thickness  is  small  relative  to  its  radius. 

The  n  =  1  mode  is  primarily  a  translational  one  which  accounts  for 
the  rigid  body  translation  of  the  shell  after  it  has  been  enveloped  by  the 
incident  wave.   Thus,  it  can  be  seen  that  the  velocity  components  for  this 
mode  approach  constant  values  equal  to  the  velocity  of  the  medium  behind  the 
incident  wave  front,  and  displacements  grow  without  bound  reaching  a  straight 
line  variation  with  time.   Note  that  the  stresses  contributed  by  this  mode 
reach  their  peak  values  within  one  transit  time  and  quickly  damp  out, 
approaching  zero  asymptotically.   For  the  incident  shear  wave,  the  n  =  0  mode 
is  also  a  rigid  body  movement  which  accounts  for  rigid  body  rotation,  and 
contributes  little  to  the  stresses. 

Modal  quantities  obtained  are  coefficients  of  Fourier  series j 
therefore,  the  total  response  or  effect  is  determined  by  adding  the  coefficients 
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multiplied  by  the  appropriate  sine  and  cosine  terms  for  any  desired  angle. 
Figures  21  through  2k   show  the  time  variation  of  stresses  in  the  shell  and 
medium  for  various  angles  when  the  first  three  modes  are  summedo   The  maxi- 
mum stress  in  the  shell  for  any  angle  may  be  determined  by  adding  the  bending 
stress  to  the  hoop  stress.   This  is  indicated  in  Figs.  21  and  2k   by  the 
dotted  line  above  the  hoop  stress. 

Figure  25  is  given  to  illustrate  the  relative  magnitudes  of  the 
hoop  stresses  in  the  shell  and  medium  for  several  thicknesses  of  shell.   This 
also  shows  the  effect  of  varying  the  relative  thickness  of  the  shell  on  the 
hoop  stress  in  the  medium »   The  dotted  line  indicating  the  hoop  stress  in 
the  medium  for  an  unlined  cavity  was  obtained  from  the  report  by  Paul  (6) . 

Figures  26  and  27  show  how  the  relative  thickness  of  the  shell 
affects  the  radial  and  shear  stresses  in  the  medium. 

As  was  discussed  earlier^  stresses  in  the  medium  for  any  radius 
can  be  determined  by  reevaluating  the  integral  terms  which  represent  the 
effects  of  the  outgoing  waves j,  and  adding  them  to  the  Fourier  expansion  of 
the  incident  wave.   Figure  28  shows  the  modal  and  total  radial,  hoop,  and 
shear  stresses  which  were  computed  for  a  time  equal  to  10  transit  times. 
These  values  are  compared  later  with  the  static  stresses,  but  on  this  figure 
the  static  stresses  do  not  differ  by  more  than  the  thickness  of  the  lines, 
and  therefore  are  not  shown c   The  time  variation  of  the  radial  and  hoop 
stresses  in  the  medium  for  various  radii  are  shown  in  Fig.  29  for  the  incident 
dilatational  wave,  and  in  Fig.  30  for  the  incident  shear  wave. 

k.J)     Short-Time  and  Asymptotic  Comparisons 

A  method  for  obtaining  a  solution  to  the  problem  which  is  accurate 
for  very  short  times  (t  «  l)  was  presented  in  Section  5»T°   This  was  desirable 
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to  validate  the  machine  solution  and  to  determine  the  effect  of  the  size  of 
time  interval  selected.   The  results  for  a  representative  problem  are  shown 
in  Figs.  51  through  53  for  the  incident  dilatational  wave,  and  in  Figs.  3^ 
through  36  for  the  incident  shear  wave.   Four  terms  of  the  series  representing 
the  F  and  G  functions,  and  three  terms  for  other  quantities  were  used  in  the 
short -time  solution. 

As  can  be  seen  from  the  graphs,  good  agreement  was  obtained  for 
very  small  values  of  time,  somewhat  shorter  time  being  obtained  for  the  shear 
wave  as  compared  to  the  dilatational  wave.   The  shorter  time  results  from  the 
nature  of  the  forcing  functions  (Eqs.  2-19  ^^^   2-32)  which  indicate  a  more 
rapid  rise  in  the  incoming  stresses  for  the  incident  shear  wave. 

Within  the  range  of  time  for  which  the  short-time  solution  is  valid, 
decreasing  the  size  of  time  interval  for  each  step  of  the  machine  solution 
resulted  in  closer  agreement  between  the  two  methods,  as  is  to  be  expected. 
It  also  indicated  that  the  stresses  and  displacements  are  not  as  sensitive  to 
variations  in  the  interval  size  as  are  the  F  and  G  functions. 

At  the  other  end  of  the  time  scale,  i.e.,  at  a  relatively  long  time 
after  passage  of  the  incident  wave  front  across  the  cavity,  another  check  on 
the  accuracy  of  the  machine  solution  is  afforded  by  the  asymptotic  approach 
of  all  values  to  the  static  results.   Figures  presented  thus  far  have  shown 
that  the  static  condition  is  approached  well  within  ten  transit  times. 

The  time  interval  used  in  the  machine  solution  affects  the  stability 
of  the  results  for  long  times.   This  is  indicated  in  Fig.  37  which  shows  the 
variations  in  computed  values  of  the  displacement  components  for  mode  n  =  2 
at  relatively  long  times,  for  different  time  intervals.  N  represents  the 
n-umber  of  time  steps  required  for  the  incident  wave  to  travel  a  distance  equal 
to  the  radius  of  the  opening,  and  1/n  defines  the  interval  size.  Note  the 
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smaller  graph  which  shows  the  percentage  difference  between  the  computed 
values  at  ten  transit  times  and  the  static  values. 

As  the  interval  is  decreased  the  machine  solution  at  long  time 
approaches  the  static  solution  more  closely.   Below  a  certain  size  of  time 
interval,  there  is  little  difference  in  the  results,  which  indicates  asymp- 
totic convergence  to  the  correct  solution.   For  this  particular  problem, 
N  -  50  seems  to  he  "critical"  in  that  for  N  <  ^0,   wide  variations  in  computed 
values  occur.  As  the  mode  number  increases,  the  "critical"  value  of  N  in- 
creases rapidly,  and  the  requirements  of  computer  storage  and  calculation 
time  become  decisive  factors  which  make  impractical  the  study  for  long  times 
of  modes  much  larger  than  2.   The  interval  size  selected  for  all  problems 
solved  (exclusive  of  the  study  to  determine  the  effect  of  the  interval  size) 
was  set  equal  to  l/i+0  (N  =  i+O)  of  the  half  transit  time  of  the  incident  wave. 

For  the  static  case,  only  the  modes  n  =  0  and  2  yield  values  other 
than  zero.   Table  1  compares  values  of  shell  stresses  and  displacement  com- 
ponents in  these  modes  obtained  from  the  computer  solution  to  a  particular 
problem  at  a  time  equal  to  10  transit  times,  with  the  static  solution.   Most 
pairs  of  values  differ  by  less  than  one  percent.   Comparable  agreement  of 
stresses  in  the  medium  at  various  radii  are  shown  in  Table  2.   Values  in 
Tables  1  and  2  were  obtained  from  the  solution  to  a  problem  whose  parameters 

were:   ti„  =  1|.0,  tj   =  3.0,  rj  =  .05,  and  v  =  .25. 
rj  p  "C 

Although  the  figures  and  tables  presented  above  were  for  a  particular 
problem,  the  discussion  given  is  applicable  to  all  problems  which  were  solved. 
Changing  the  physical  characteristics  of  the  shell  and  medium  within  the  range 
of  values  studied  had  hardly  any  effect  on  the  degree  to  which  the  long  time 
machine  solution  and  the  static  solution  agreed.  Also,  the  absolute  value  of 
all  quantities  which  should  asymptotically  approach  zero  became  less  than 
.00005  in  each  case  well  before  ten  transit  times. 
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k.k-     Effect  of  Parameters 

Studies  were  conducted  to  determine  the  effect  of  each  of  the 
following  parameters : 


T]  =  —        V   =  V  (medium) 


for  an  unstiffened  shell  without  additional  mass.   Since  it  was  impractical 
to  take  into  account  all  permutations  of  the  above  parameters,  a  basic  shell 
where 

n^  =  .05  V    =  .25 

was  considered  from  which  each  parameter  was  separately  varied  to  determine 
its  effect  on  the  resulting  stresses  and  displacements.   Calculations 
performed  were  only  for  the  case  of  the  incident  dilatational  (P)  wave. 

Of  particular  interest  was  the  determination  of  the  maximum  dynamic 
stresses  and  displacements  (not  including  rigid  body  translation)  due  to  the 
incident  stress  wave.   Tables  3  through  6  compare  the  maximum  values  obtained 
in  the  machine  solution  with  the  static  solution.   "DLF/'  termed  the  dynamic 
load  factor,  is  defined  as  the  factor  by  which  the  displacement  or  stress 
produced  by  a     applied  as  a  static  load  should  be  multiplied  in  order  to  obtain 
the  maximum  dynamic  value.   Figures  38  through  ^1  are  given  to  graphically 
illustrate  the  variations  in  stresses  and  displacements,  both  static  and 
dynamic,  in  the  range  of  parameters  considered.   Stresses  are  given  in  units 
of  Ja  I  and  displacements  in  non-dimensionalized  units  defined  by  equations 
(2-36)  and  (2-50). 
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Table  3  and  Fig.  38  illustrate  the  effect  of  increasing  the 
relative  thickness  of  the  shell..   The  range  of  values  selected  for  r\      is 
probably  much  greater  than  is  practical  or  justified  by  the  assumptions  of 
the  analysis,  and  was  considered  only  to  determine  the  trend  of  the  results. 
As  T]   is  increased,  displacements  and  hoop  stresses  of  both  the  shell  and 
medium  decrease;  however,  the  bending  stress  of  the  shell  and  the  radial  and 
shear  stresses  in  the  medium  at  the  boundary  increase.   Figure  25  shows  the 
time  variation  of  the  hoop  stress  in  the  shell  and  medium  for  various  thick- 
nesses of  shell,  including  the  case  of  the  unlined  cavity. 

Increasing  t)  results  in  a  rapid  increase  in  the  displacements  and 
stresses  of  the  shell  and  a  much  lower  rate  of  attenuation  of  the  hoop  stress 
in  the  medium.   See  Fig.  39  a-^^d.  Table  k. 

Table  5  and  Fig.  ^0  show  that  the  parameter  t]  has  no  effect  on  the 
static  results  but  does  affect  the  dynamic  response.  Increasing  r\  increased 
the  maximum  response  in  both  the  medium  and  shell. 

Figure  ^1  and  Table  6  show  that  within  the  range  of  Poisson's  ratio 
for  the  medium  considered,  as  v  increases  displacements  and  stresses  in  the 
shell  decrease,  with  very  little  additional  reduction  of  the  medium  hoop 
stress . 

Additional  mass  within  the  shell  was  shown  to  contribute  relatively 
little  to  the  overall  response.   This  is  mainly  due  to  the  simplifying  assump- 
tions which  were  made  in  deriving  the  equations  of  motion  to  account  for  the 
additional  mass.   Figure  k2   shows  the  effect  of  the  presence  of  additional 
mass  equal  to  k-O-Tr   times  the  mass  of  the  shell  on  the  dynamic  response  of  a 
particular  shell.  Additional  mass  decreases  the  displacements  and  increases 
the  dynamic  stresses  in  mode  n  =  1. 
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ifo5  Response  to  Time  Varying  Incident  Wave 

The  results  obtained  from  a  solution  to  an  incident  wave  with  a  step 
variation  in  time  was  shown  to  be  useful  by  application  of  Duhamel's  integral  in 
determining  the  response  to  any  time  varying  stress  wave.   The  case  of  the 
exponentially  decaying  wave  was  considered.,  the  results  of  which  are 
illustrated  in  Fig.  43-   Substantial  reductions  in  the  maximum  stresses  can 
be  expected  as  the  duration  becomes  smaller.   In  Fig.  kk-   is  shown  the  effect 
of  a  linear  rise  in  the  amplitude  of  the  incident  wave.   Note  that  for  a 
wave  with  a  linear  rise  followed  by  a  step  variation  in  time,  very  little 
decrease  in  maximum  stresses  occurred. 

k.6     Comparison  with  Previous  Work 

Baron  (l),  using  a  different  method  of  analysis,  investigated  the 
dynamic  response  of  two  shells  with  different  physical  characteristics  sub- 
jected to  an  incident  plane  dilatational  wave.   Figures  ^5  through  ^9  compare 
his  results  for  the  modal  values  of  stresses  and  displacements  of  the  thin 
shell,  with  results  obtained  by  the  method  of  solution  given  in  this  study » 
Similar  results  were  obtained  considering  his  so-called  stiff  shell..  Although 
the  shape  of  the  response  curves  can  be  considered  similar  in  the  two  reports, 
the  magnitude  of  the  dynamic  response  in  Ref.  (l)  seems  consistently  higher 
than  in  the  present  report,  and  the  long  time  results  asymptotically  approach 
values  higher  than  the  static  solution.   The  following  table  compares  the 
maximum  modal  stresses,  for  Poisson's  ratio  of  0 .25> obtained  in  the  two  reports, 
and  also  shows  values  for  the  static  case.   Stresses  are  given  in  units  of 
I  a  j .   The  bending  stresses  for  the  stiffened  shell  cannot  be  compared  directly 
since  d,  the  distance  from  the  neutral  axis  of  the  shell  to  its  extreme  fiber, 
is  not  stated  in  Refo  (l) . 
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Thin  Shell 

Stiff  Shell 

Maximum 

stresses 

Maximum 

Stresses 

Ref. 

Current 

Ref. 

CiJ.rrent 

Quant . 

Mode 

(1) 

Work 

Static 

(1) 

Work 

Static 

^99 

0 

-1+.T9 

-i+.56 

-I^-.IO 

-4.79 

-i|.08 

-3. 67 

^99 

2 

h.G9 

1^.11 

3.67 

3A2 

3.20 

2.95 

"^sb 

0 

o.oi^3 

0.01^3 

0.039 

0.55 

l|.08d 

3.67d 

^sb 

2 

O.li^O 

0.137 

0.122 

1,70 

13 . 29d 

11 . 89d 

^^-•7  Conclusions 

Conclusions  drawn  from  the  results  of  the  analysis  are: 

(1)  The  method  which  has  "been  presented  is  practical  for  effective 
computation  of  the  dynamic  response  of  a  cylindrical  shell  embedded  in  an 
elastic  medium  when  subjected  to  plane  dilatational  or  shear  waves.   The 
solution  presented  herein  is  believed  correct  since  it  was  checked  by  independent 
methods  at  short  and  long  times. 

(2)  Peak  stresses  and  displacements  in  both  the  medium  and  shell 
occur  sometime  after  the  transit  of  the  incident  wave  across  the  cavity;  of 
the  problems  solved,  both  the  average  and  mean  time  at  which  the  peak  values 
occurred  was  equal  to  3  transit  times.   The  dynamic  effect  measured  in  terms 
of  the  ratio  of  the  maximum  stresses  and  displacements  to  the  static  values 
varied  within  a  relatively  small  range.   The  average  value  of  this  ratio  was 
I0I2  and  the  mean,  1.11.   Thus,  for  the  cases  considered  at  least,  the  maximum 
stresses  and  displacements  to  be  expected  for  a  particular  situation  can  be 
roughly  approximated  by  determining  the  static  values  and  by  multiplying  them 
by  a  factor  of  1.1. 
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(5)  The  largest  stresses  in  the  shell  occurred  for  a  relatively 
thin  liner  in  a  medium  with  a  low  modulus  of  elasticity,  low  mass  density, 
and  low  Poisson's  ratio.   The  greatest  reduction  of  the  hoop  stress  in  the 
medium  as  compared  to  the  unlined  cavity  results  from  a  relatively  thick 
liner  in  a  medium  with  low  modulus  of  elasticity  and  high  mass  density. 
Additional  mass  within  the  shell  has  relatively  small  effect  on  the  dynamic 
stresses . 

(k)     The  practical  value  of  tunnel  linings  to  reduce  the  maximum 
stresses  in  the  medium  depends  on  the  several  conditions  mentioned  above,  and 
on  the  magnitude  of  the  incident  stress  wave.   Under  favorable  conditions, 
reduction  of  stresses  on  the  order  of  30  percent  or  more  is  possible.   However, 
for  materials  such  as  granite,  smaller  reductions  can  be  expected  for  steel 
liners  of  practical  size. 

(5)   Certain  assumptions  made  concerning  the  behavior  of  the  shell 
have  limited  the  applicability  of  the  analysis  to  relatively  thin  liners. 
Future  studies  of  the  behavior  of  thick  shells  would  indicate  the  effect  of 
the  approximations  used  herein.   The  analysis  has  also  been  based  on  the 
assumption  that  the  behavior  of  both  the  liner  and  medium  is  linearly  elastic 
at  all  times.   Perhaps  a  more  desirable  condition  would  be  one  in  which  some 
inelastic  behavior  is  permitted  to  take  place  in  the  medium  surrounding  the 
shell  or  in  the  shell,  or  one  in  which  some  inelastic  energy  absorbing  medium 
such  as  cinders  or  foamed  plastic  surrounds  the  shell.  As  a  subject  of  future 
study,  it  is  recommended  that  the  behavior  of  thick  shells,  and  thin  shells 
surrounded  by  some  energy  absorbing  layer,  be  considered. 
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APPENDIX  A 
DERIVATION  OF  THE  EXPRESSIONS  INVOLVING  THE  POTENTIAL  FUNCTIONS 

A.l  General  Form  of  the  Potential  Function 

The  dilatational  potential  function  must  satisfy  the  wave 
equation 

^+i^  +  ^^-i_^  (A-i) 

A  2   r  ^r  ^  2  .^2  -  2  .^2  ^^  ^^ 

OT  T       d^  C   Ot 

cp  in  the  case  of  the  incoming  dilatational  wave  is  an  even  function  of  9   and 
thus  can  be  expressed  as  a  cosine  Fourier  series 

QQ 

n=0 
where  f  is  a  function  of  r  and  t  representing  the  modal  coefficient  of  the 
potential  function.   By  substitution  of  the  above  expression  into  equation 
(A-l)  we  obtain  the  equation  that  must  be  satisfied  by  f  . 

a^f      ,  af       2         ,  a^f 

.  2  ^  r  ^r     2  n  ~  2  .,2  ^^  -'^ 

c3r  r       c,  ot 

The  general  solution  of  f  is  assumed  to  be  of  the  form 
to  n 

f  =  r"^  R  (r,t)  (A-k) 

n       n^  -^  '  \    / 

Then  R     is   some   function  of  r  and  t  satisfying 


2^  -s2 
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at^       1 


a  R          /„     -,  N   ^R  " 
n        (2n+lj   n 

2  r         ^r 


-b2 

If  R  is  a  solution  it  can  be  shown  that  the  corresponding  equation  for 


(A-5) 


n 

R  -,  is  satisfied  by 
n+1  "^ 


^R 
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By  repeated  application  of  this  result  it  can  be  shown  that  equation  (A-5) 
is  satisfied  by 

R  .  (i  !_)''  R  (A-T) 

n    ^r  or    o 

where  R  is  the  solution  of 
o 

.2„        ^2. 
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rd  R  ,  Sr 

-^  +  -^1                       (A-8) 

-.  2  r  dr  ' 
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The  solution  to  this  equation  for  the  case  of  a  wave  diverging  from  a  center 


IS 

00 


R  =   /  '  F(t  -  —  cosh  u^  )  du  (A-9) 

o    Jq       c^       -    1 

For  the  proof  of  this,  see  Lamb  {h) .      Therefore,  the  coefficient  f  of  the 
potential  function  is  written  as 

op 

f  =  r'^(i  |_)    /•  F(t  _  i—  cosh  uj  du,  (A-IO) 

n     ^r  dr    J  c  1'    1 

For  modes  0,  1,  and  2,    this  can  be  written 


00 

00 


•1^0 

00  00 

fr,  =^     I     F"(^)  coshV  du  +;--/■  F'(f  )  cosh  u  du 
"^c^^O  ^^^'^l^O 

where  primes  indicate  the  derivative  of  the  function  F  with  respect  to  its 
argument.   Note  that  the  expression  for  the  n  =  2  mode  contains  derivatives 

of  two  different  orders.   It  is  convenient  to  express  the  function  f  in 

^  n 

terms  of  derivatives  of  a  single  order  which  may  be  done  through  integration 
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"by  parts.  For  the  n  =  2  mode  we  then  wish  to  change  the  second  integral  to 
an  integral  involving  the  second  derivative  of  the  function  F.  Integrating 
by  parts  gives 


u 


Im 


^dT)   =   i'T] 


^Im 


u 


Im 


•q-d^ 


where 


Thus 


i  =  F-a^) 


d|  =  -  f-  F"(li)  sinh  u^  du^ 


T]  =  sinh  u 


1 


dri  =  cosh  u,  du. 


00  00  U 

/■    F'(^    )    cosh  u     du     =—      /■    F"(|    )    sinh^u  du     +     F' (^    )    sinh  u 
Jq  1  IIc^Jq  1  11        _l  IJ^ 


The  second  term  can  he  shown  to  equal  zero  since  for  its  lower  limit 
sinh  U-,  =  0;  and  for  its  upper  limit,  F'  is  the  integral  of  the  function  F 
at  the  wave  front  which  under  the  assumed  initial  conditions^  does  not  exist. 
Therefore,  the  coefficient  of  the  potential  function  for  the  n  =  2  mode  can 
be  written  as 


1 


f  =  — -   /   F"(t  -  — ■  cosh  u  j  cosh  2u  du 


(A -11) 
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By  a  similar  process  the  coefficients  of  the  potential  function  for  any  other 
mode  can  be  reduced  to  the  form 


ill) 
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00 

-^  /   F  (t  -  — —  cosh  U-,  )  cosh  nu^  du- 

n   J  ^    c        1'  1        . 

1 


(A-12) 


The  shear  potential  function  must  satisfy  the  wave  equation 

b%    1  at    1   ^%     1  d% 


ar^  ■  ^  ^  ■  r^  ^^ 


2  ;^+2 
Cg  dt 


(A -13) 
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i|f  in  the  case  of  the  incoming  dilatational  wave  must  "be  an  odd  function  of  9 
and  can  te  represented  as 


lif  =  )   g^(r,t)  sin  ne 


n=l 
and  proceeding  exactly  as  in  the  case  of  the  dilatational  potential  the 
general  expression  for  the  coefficient  g  is  found  to  be 

00 

g  =  i^ii-  /  G^{t    -  —   cosh  u^)  cosh  nuo  du^  (A-lU) 

n     ^n   Jq   ^    c^       2'  2   2 

A. 2  Velocity  Terms 

The  equations  for  the  velocity  components  of  a  particle  on  the 
"boundary  due  to  the  outgoing  waves  are 

:.2     ,  ^2 

""  ~  "Mt  ^  r  ae§t 

.  .2     ,2,  (^-15) 

^  '  r  d9dt  "  ^rSt 

The  components  of  velocity  and  the  potential  functions  are  expanded  in  series 


as 


00 

=     )      a    (r^t)    cos  nS 


u 

n=0 


Y  ^     )      b    (r^t)    sin  nd 

n=l 

(A-16) 


=     )      f   (r^t)    cos  xiQ 


9 

n=0 


^  =     /      g^(r,t)    sin  ne 
n=0 
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Substitution  into  equations  (A-.15J  gives 

^^f      ^g 

.  /   ,  \      n   n  °n 

P  (A-IT) 

^f    a  g 
C  /   +^     n   n      n 

The  expressions  for  the  coefficients  f  and  g  are  given  in  equations  (A-12) 
and  (A-1^-).   Using  these  in  the  above  equations  and  by  the  application  of 
integration  by  parts  we  get  the  following  velocity  terms 

CO 
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A. 3  Stress  Terms 

Stresses  in  the  medium  due  to  the  outgoing  waves  are  written  in 
series  form  as 

00 


o        =  )  A  (r.t)  cos  n0 
n=0 
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By   substitution  into  equations  (2~10j  the  coefficients  A  ,  B  ,  and  C  are 
expressed  in  terms  of  the  coefficients  of  the  potential  functions,  f  and 
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B  (out)  =  u. 
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Again„  using  tiie  equations  for  f  and  g  given  earlier  we  can  write  the  above 
coefficients  finally  as 
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APPENDIX  B 
STATIC  SOLUTION 

B.l  Dilatational  Wave 

The  static  solution  presented  here  is  based  on  the  application  of 
the  theory  of  stress  functions  presented  by  Timoshenko  (9) •   Under  static 
conditions^  it  can  be  assumed  that  at  large  distances  from  the  boundary  of 
the  cavity  the  state  of  stress  in  the  medium  is  equal  to  the  stress  field  in 
the  medium  behind  the  front.   In  polar  coordinates  this  is 

a   =  -  (i±^)  a  -  (i^)  a     cos  29 
rr       2    p     2    p 

"re   =  (^)  °p  =^"  29  (B-1) 

"ee  '   (¥>  "p  ^  (¥)  °p  -=  29 

a   is  the  stress  in  the  medium  in  the  direction  of  wave  propagation  and  vo 
P  ^  ^  ^  p 

is  the  stress  parallel  to  the  wave  front,   v^  derived  from  the  assumption 
that  there  is  no  strain  parallel  to  the  wave  front,  is  equal  to 

V   =  T^  (B-2) 

1-v 

It  can  be  seen  from  equations  (B-1)  that  the  n  =  0  and  2  modes  describe 
exactly  the  free  field  stresses.   Therefore,  the  unknown  stresses  at  the  edge 
of  the  cavity  can  likewise  be  expressed  in  terms  of  these  modes 


a   =  -  P  -  P^  cos  29 
rr      o    2 


\q   =  Sg  sin  29 


(B-5) 


where  P  ,    Pp^  and  S  are  the  unknown  modal  components  of  stress  acting  at 
the  boundary.   These  same  stresses  must  act  on  the  shell.   Thus  the  boundary 
stresses  in  the  medium  and  shell  can  be  illustrated  as  shown  in  Fig,  9* 
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The  equations  of  equilibrium  in  the  medium  under  plane  strain 
conditions  are  satisfied  by  the  following  expressions  for  the  components  of 
stress 

_  i  M   A-  ^ 
rr  ~  r  ^r    2  >,„2 
r  ad 


1  d(2   1  6^ 


t9  2  bd        r  brb9 
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_  ^ 

^ed      .  2 

dr 


(B-i^) 


where  ft  is  the  stress  function  in  terms  of  r  and  Q. 

B.ll  n  =  0  Mode 

The  general  solution  of  the  stress  function  for  the  n  =  0  mode  is 
given  by  Timoshenko  as 

2 
fi  =  K  log  r  +  Mr 

from  which  can  be  derived  the  stress  components 
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The  coefficients  K  and  M  are  determined  from  the  states  of  stress  at  r  — ^  °o 
and  r  =  R 


K  =  R 


<¥)  "p  -  ^0 


^       k      ^  P 

For  the  conditions  of  plane  strain^  the  strains  are 


(B-6) 
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where 
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modulus  of  plane  strain 
for  the  medium 


Displacement  components  are  found  "by  suitable  integration  of  the  following 
equations  for  the  strains  in  the  medium 
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Displacements  in  the  medium  for  the  zero  mode  are 
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The  corresponding  displacement  components  of  the  shell  with  an  exterior 
compressive  force  P  are  found  in  Flugge  (2)  to  be 
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where 
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modulus  of  plane  strain 
for  the  shell 
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The  unknown  stress  P  is  now  determined  by  equating  the  radial  displacement 
components ;  to  get 


(1+v)  o 
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(B-11) 


After  P  is  determined^  stresses  in  the  medium  at  any  radius  can  "be  found 
using  equations  (3-5)"   Equation  (B-10)  gives  the  displacement  of  the  shell; 
and  the  hoop  and  bending  stresses  of  the  shell  are^  respectively 
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B.12  n  =  2  Mode 


The  general  solution  in  terms  of  a  stress  function  Q   is  given  as 
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from  which  the  stress  components  in  the  medium  become 


rr 


2K  +  -j^  +  — 
r    r  -■ 


cos  29 


r  -' 


re 


2K 


cos  29 
6M   2N 


r    r  -• 


sin  29 


(B-13) 


Coefficients  are  determined  as  before  from  the  states  of  stress  at  r  — >  <» 
and  r  =  R 


K  =  {^)   a 
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(B-11+) 
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M     = 

;2%^  "p  -  p,  -  SB,; 

N      = 

r2 

2 

_-(l-v)   a^  .  P^  +  S2_ 

The  displacement  components  of  the  medium  at  the  boundary  can  be  expressed 


as 
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and  the  corresponding  components  for  the  shell  are 
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The  unknown  boundary  stresses  P  and  S  are  determined  by  equating  the  dis- 
placement of  the  shell  and  medium  at  the  boundary.   The  hoop  and  bending 
stresses  in  the  shell  are  then  found  from  the  following  equations 


dQ 


SB 


PR   2S  R 
(_  +  __)  eos  29 
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j^   (2P2  +  Sg)  cos  29 


(B-17) 


Displacements  of  the  shell  and  stresses  in  the  medium  can  be  determined  usxng 
equations  given  earlier  in  the  discussion.   Total  static  stresses  and  dis- 
placements are  merely  the  sum  of  the  modal  values  for  any  angle. 
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>.2  Shear  Wave 


The  free  field  shear  wave  stresses  expressed  in  polar  coordinates 
include  only  the  n  =  2  mode 


a 

rr 

a     sin  20 
s 

"tQ  = 

a     cos   29 
s 

""qq  = 

-a     sin  29 
s 

(B-18) 


where  a  is  the  amplitude  of  the  incoming  wave.   The  static  solution  is 
obtained  exactly  as  in  the  case  of  the  dilatational  wave  except  for  the 
interchange  of  sines  and  cosines  resulting  from  the  difference  in  geometry. 
The  resulting  equations  are 

a.   Displacement  components 
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"b.  Shell  Stresses 
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11^  (2P2  +  S^)  sin  20 


(B-20) 


c.  Medium  Stresses 


a   =  _  (2A  +  ^  +  it§)  sin  20 
rr     ^      1+2' 
r    r 


r0 


00 


r    r 

(2A  +  ^)  sin  20 
r 


(B-21) 


where 
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A  =  -  - 
.2 


D 


R 


2°s  -  ^2  -  ^2 


In  the  machine  program,  the  n  =  0  mode  is  purely  rotational  and 
thus  we  get  a  static  value  of  tangential  displacement.   This  is  calculated 
as  follows 


9    =i7    =ii 
rot   2  xy   2  |j. 


(B-22) 


5  =  « 


1  's 


rot   2  |j. 
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DILATATIOHAL  WAVE 


time  =  20  R/ci 


Node 

Quant. 

Stat.     Mach. 

Dev.   From 
Static 

Mode 

Quant 

Stat. 

Mach. 

Dev.   Prom 
Static 

0 

-oL 

4.000    »>.oo6 

.15^ 

2 

-OC 

4.613 

4.599 

.30^ 

0 

-Oie 

if.OOO     U.006 

.15 

2 

/3 

3.845 

3.856 

.29 

0 

(Tsb 

.100        .100 

.00 

2 

CToe 

3.077 

3.114 

1.20 

2 

<r»b 

.3^6 

.3»*5 

.29 

Stresses  gl>ren  in  units  of  \<rp\ 

Oisplaceaents  defined  by  Bqs.  (2>50)  and  (2-36) 


SHEAR  WAVE 


time  =  20  R/c2 


Dev.  FroBi 

Node 

Quant. 

Stat.     Mach. 

Static 

0 

(S 

5.000     4.996 

.08^ 

2 

OL 

13.840  13.930 

.65 

2 

/3 

11.535  11.539 

.03 

2 

-CJbb 

9.232     9.1^9 

.90 

2 

CTsb 

1.038     1.045 

.67 

TABLE  1 
CONPARISOH  or  STATIC  AND  MACHINE  VALUES  FOR  MODAL  RESPONSE  OF  THE  SHELL 
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DILATATIOHAL  WAVE 


time  =  20  R/ci 


Mode 

Quant* 

1.0 

R 

2.0 

R 

3.0 

R 

4.0 

R 

Stat. 

Mach. 

Stat. 

Mach. 

Stat. 

Mach. 

Stat. 

Mach. 

0 

-  (TfY 

.200 

.200 

.550 

.551 

.615 

.616 

.638 

.639 

0 

-CTee 

1.133 

1.135 

.783 

.785 

.718 

.720 

.696 

.697 

2 

-(Trr 

-.153 

-.155 

.111 

.112 

.226 

.229 

.271 

.273 

2 

OVe 

.307 

.312 

.428 

.429 

.384 

.386 

.363 

.366 

2 

OeB 

.872 

.868 

.367 

.368 

.340 

.342 

.335 

.337 

*  Stresses  given  in  units  of  lfp| 


SHEAR  WAVE 

time  =  20  R/c2 

Mode 

Quant. 

1.0  R 

2.0  R 

3.0 

R 

4.0  R 

Stat.  Mach. 

Stat.  Mach. 

Stat. 

Mach. 

Stat.  Mach. 

2 

(Ty^ 

-.460  -.455 

.332   .331 

.675 

.678 

.814   .815 

2 

CJVe 

.921   .913 

1.283  1.286 

1.151 

1.154 

1.090  1.092 

2 

-Oee 

2.616  2.636 

1.101  1.104 

1.020 

1.022 

1.006  1.009 

TABLE  2 
CQMPARISOH  OF  STATIC  AND  MACBIRE  VALUES  PQR  MODAL  STRESSES  IH  THE  MEDHI4 
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\ 

^Be,  (90-) 

DIZ* 

t„# 

^^^ 

(90-) 

mr 

tB 

static 

Max. 

Static 

Max. 

.01 

-9.21         -10.12 

1.10 

3.12 

-2.'*9 

-2.73 

1.09 

3.13 

.05 

-7.08         -  7.7^ 

1.09 

2.88 

-2.01 

-2.19 

1.09 

2.88 

.10 

-5.56        -  6.15 

1.10 

2.25 

-1.67 

-1.84 

1.11 

2.25 

.15 

-U.60        -  5.21 

1.13 

2.25 

-1.45 

-1.65 

1.13 

2.25 

.20 

-3.9**         -  h.52 

l.llf 

2.25 

-1.30 

-1.50 

1.15 

2.25 

\ 

C»b 

(90") 

mr 

ta 

Crvr*y> 

(90-) 

ULF 

tm 

Static 

Max. 

Static 

Max 

.01 

.10 

.11 

1.11 

3.38 

-  .09 

-    .10 

1.10 

3.25 

.05 

M 

.50 

1.11 

3.25 

-  .35 

-  .39 

1.11 

2.88 

.10 

.83 

.9** 

1.13 

3.00 

-  .55 

-  .63 

l.llf 

2.50 

.15 

1.18 

1.36 

1.15 

2.88 

-  .68 

-  .79 

1.17 

2.50 

.20 

1.50 

1.76 

1.18 

2.88 

-  .75 

-  .91 

1.21 

2.50 

\ 

OL 

(0-) 

DLF 

tm 

/3 

(45') 

DLF 

ta 

Static 

Max. 

Static 

Max. 

.01 

-9.65 

-10.64 

1.10 

3.12 

4.67 

5.23 

1.12 

3.75 

.05 

-8.61 

-  9.50 

1.10 

3.00 

3.85 

4.33 

1.13 

3.50 

.10 

-7.76 

-  8.62 

1.11 

2.76 

3.33 

3.82 

1.15 

3.25 

.15 

-7.15 

-  8.01 

1.12 

2.88 

3.02 

3.55 

1.18 

3.25 

.20 

-6.66 

-  7.56 

1.13 

2.88 

2.80 

3.37 

1.20 

3.38 

*  Ratio  of  Max.  to  Static  #  Time  in  transit  times  at  vhich  Max.  occurs 

TABLE  3 
EFFECT  OF  RELATIVE  THICKMESS  OF  SHELL 
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\ 

Osftj 

(90*) 

Diy* 

t/ 

<J»e>. 

(90*) 

— 1 

DLF 



tn 

Static 

Max. 

Static 

Max. 

2 

-  U.16 

-  ^.57 

1.10 

2.76 

-2.27 

-2.51 

1.11 

2.76 

k 

-  7.08 

-  7.7^ 

1.09 

2.88 

-2.01 

-2.19 

1.09 

2.88 

6 

-  9.33 

-10.11 

1.08 

2.76 

-1.81 

-1.97 

1.08 

2.76 

10 

-12.57 

-13.1*1 

1.07 

2.76 

-1.55 

-1.65 

1.07 

2.76 

16 

-15.72 

-16.47 

1.05 

2.88 

-1.31 

-1.37 

1.05 

2.76 

r\ 

OVb 

(00**) 

<^yy^ 

fgo") 

\ 

static 

Max. 

DIF 

ta 

static 

Max. 

DU 

tm 

2 

.23 

.26 

1.12 

3.25 

-  .21 

-  .23 

1.13 

3.12 

k 

.**5 

.50 

1.11 

3.25 

-  .35 

-  .39 

1.11 

2.88 

6 

.6k 

.71 

1.11 

3.12 

-  M 

-  .51 

1.10 

2.76 

10 

1.01 

1.12 

1.10 

3.12 

-  .63 

-  .68 

1.08 

2.62 

16 

1.53 

1.68 

1.09 

3.12 

-  .78 

-  .83 

1.06 

2.62 

OL 

(0-) 

mF 

tm 

(S 

(45*) 

DI^ 

ta 

Static 

Max. 

Static 

Max. 

2 

-  4.60 

-  5.12 

i.n 

3.00 

2.14 

2.43 

1.13 

3.62 

k 

-  8.61 

-  9.50 

1.10 

3.00 

3.85 

4.33 

1.13 

3.50 

6 

-12.22 

-13.37 

1.09 

2.88 

5.32 

5.96 

X.12 

3.38 

ID 

-18.65 

-20.20 

1.08 

2.88 

7.92 

8.82 

1.11 

3.25 

16 

-27.11 

-29.15 

1.08 

3.00 

11.41 

12.64 

1.11 

3.12 

'Batio  of  Max.  to  Static 


'Tine  in  transit  tlaes  at  which  Max.  occurs 
TABLE  4 
KPPECT  OF   RATIO  OF  MODULI  OF   ELASncmr 
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> 

o"w* 

(90^) 

DLF* 

t»# 

Oie^ 

(90*) 

Dur 

tn 

Static 

Max. 

static 

Max. 

1 

-7.08 

-7.67 

1.08 

3.00 

-2.01 

-2.17 

1.08 

2.88 

3 

II 

-7.7^ 

1.09 

2.88 

(1 

-2.19 

1.09 

2.88 

5 

II 

-7.97 

1.13 

2.50 

11 

-2.26 

1.12 

2.50 

7 

It 

-8.30 

1.17 

2.61 

II 

-2.3*^ 

1.17 

2.62 

10 

n 

-8.70 

1.23 

2.88 

II 

-2.1f5 

1.22 

2.88 

\ 

a$b 

(90*) 

DLT 

tm 

Crrw. 

i^") 

mr 

1 
tm 

static 

Max. 

Static 

Max. 

1 

M 

.^9 

1.10 

3.25 

-  .35 

-  .38 

1.09 

3.00 

3 

H 

.50 

1.11 

3.25 

M 

-  .39 

1.11 

2.88 

5 

n 

•  50 

1.13 

3.12 

II 

-  .ki 

1.15 

2.88 

7 

n 

.51 

1.15 

3.25 

11 

-  .^3 

1.21 

3.12 

10 

II 

.53 

1.18 

3.50 

II 

-  M 

1.30 

3.50 

> 

CL 

(0^ 

DUr 

tn 

/3 

(k5^) 

DLF 

tn 

Static 

Max. 

Static 

Max. 

1 

-8.61 

-9.38 

1.09 

3.00 

3.85 

U.27 

1.11 

3.50 

3 

» 

-9.50 

1.10 

3.00 

11 

'♦.33 

1.13 

3.50 

5 

It 

-9.6»f 

1.12 

2.88 

II 

k,kl 

1.15 

3.62 

7 

It 

-9.77 

l.U 

3.00 

It 

kM 

1.16 

2.76 

10 

tt 

-9.95 

1.16 

2.62 

II 

4.72 

1.23 

3.00 

'Ratio  of  Max.  to  Static 


le  in  transit  times  at  which  Max.  occurs 
TABLE  5 
SFFECT  OP  MASS  RATIO 
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t) 

CJeev 

(90') 

DLF* 

ij 

0"e9>n 

(go^ 

DLF 

tm 

Static 

Max. 

Static 

Max. 

0 

-  8.33 

-  9.17 

1.10 

3.25 

-2.08 

-2.29 

1.10 

3.25 

.1 

-  7.96 

-  8.69 

1.09 

3.12 

-2.06 

-2.2U 

1.09 

3.12 

.2 

-  7.^2 

-  8.08 

1.09 

3.00 

-2.02 

-2.20 

1.09 

3.00 

.25 

-  7.08 

-  7.7V 

1.09 

2.88 

-2.01 

-2.19 

1.09 

2.88 

.k 

-  5.76 

-  6.75 

1.17 

3.12 

-1.91 

-2.22 

1.16 

3.25 

V 

(Tsb 

(90*) 

DI^ 

tm 

C5v^m 

f90') 

DLP 

t. 

Static 

Max. 

Static 

Max. 

0 

.61 

.69 

l.lU 

3.25 

-  A2 

-  .^7 

1.13 

3.00 

.1 

.56 

.63 

1.13 

3.25 

-   .ItO 

-  .kk 

1.12 

2.88 

.2 

.^9 

.55 

1.12 

3.25 

-  .37 

-  .kl 

1.11 

2.88 

.25 

.^5 

.50 

1.11 

3.25 

-  .35 

-  .39 

1.11 

2.88 

.k 

.27 

.31 

l.lU 

3.75 

-  .29 

-  .33 

1.15 

3.88 

V 

oc 

(0-) 

DIF 

ta 

fB 

(k5'') 

TSUF 

tm 

static 

Max. 

static 

Max. 

0 

-10.33 

-ll.i»8 

1.11 

3.25 

6.00 

6.9»> 

1.16 

3.38 

.1 

-  9.88 

-10.90 

1.10 

3.12 

5.31 

6.08 

1.14 

3.38 

•2 

-  9.12 

-10.03 

1.10 

3.00 

k.ko 

U.98 

1.13 

3.38 

.25 

-  8.61 

-  9.50 

1.10 

3.00 

3.85 

V.33 

1.13 

3.50 

.k 

-  6.5^ 

-  7.61 

1.16 

3.50 

1.77 

1.97 

1.11 

1^.38 

Ratio  of  Max.  to  Static 


'^TiJBe  in  transit  tines  at  which  Max.  occurs 
TABLE     6 
EFFECT  OF  POISSOH'S  RATIO 
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(Tee 


i 


0>  -[J— «f 

t 


Plane  Dllatatlonal  Wave  Front  at  time  t  -  0 
Stresses  Behind  Wave  Front  Sbovn  In  Both  Rectangular 
and  Polar  Coordinates 


Plane  Dllatatlonal  Wave  Front  at  time  t>0 
FIGURE  1 
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060 


(Trr    ^©/ 


-// 


/ 


/A 


D 


CTs 


Pleme  Shear  Wave  Front  at  time  t  =  0 
Stresses  Behind  Wave  Front  Shown  in  Both  Rectangular 
and  Polar  Coordinates 


Plane  Shear  Wave  Front  at  time  t  >  0 


FIGURE  2 
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FIG.  3  DISPLACEMENT  COMPONEIfTS  OF  THE  SHELL 


Wave  front 


FIO.  k     DEFIHITIOH  OF  POSITION  ANGLE,  Oi 
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xy  plane 


FIG.  3  GBOMETRI  OP  A  POIHT  SOURCE  ON  IBB  z  AXIS 


xy  pleme 


PIG.  6  GEOMETRY  OF  MANY  POINT  SOURCES 
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Wave  front 


FIG.  7  POSITION  ANGLE,  0,  ,  FOR  ANY  RADIUS 


S> 


*PO 


\ 

\^ 

^^v^ 

^"'^^ 

""^^.^ 

"■^^.^ 

FIG.  8  TIME  DEPENDENT  STRESS  WAVE 
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<r^ 


FIG.  9  STATIC  STRESSES  ACTING  ON  THE  MEDIUM  AND  SHELL 
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VITA 


Takeshi  Yoshihara  was  'born  Ik  May  1951  in  Selleck,  Washington ^ 
and  spent  his  early  years  in  the  Pacific  Northwest,  including  three  and 
one  half  years  internment  in  a  War  Relocation  Center.   In  June  19^9?  after 
graduation  from  high  school,  he  became  the  first  Japanese -American  to  be 
accepted  by  the  United  States  Naval  Academy  from  where  he  graduated  in  1953^ 
receiving  the  Bachelor  of  Science  degree  "with  distinction."   He  was  com- 
missioned Ensign,  United  States  Navy,  and  has  been  on  continuous  active  duty, 
presently  holding  the  rank  of  Lieutenant  in  the  Navy  Civil  Engineer  Corps. 
In  1955  he  was  sent  to  Rensselaer  Polytechnic  Institute  for  further  studies 
and  in  195^  received  the  Bachelor  of  Civil  Engineering  degree  from  that 
institution.   In  I96O  he  was  ordered  to  the  University  of  Illinois  as  a 
graduate  student  from  where  he  received  the  Master  of  Science  degree  in 
1961, 

He  is  a  member  of  the  following  honorary  societies-   Chi  Epsilon, 
Phi  Kappa  Phi,  Sigma  Xi,  and  Tau  Beta  Pi„   He  is  also  an  Associate  Member 
of  the  American  Society  of  Civil  Engineers,  and  is  a  licensed  professional 
engineer  in  the  State  of  Illinois. 
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